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Abstract 



o 



Aims. Excitation of far-infrared and submillimetric molecular lines may originate from nonreactive collisions, chemical formation, 
or far infrared, near-infrared, and optical fluorescences. As a template, we investigate the impact of each of these processes on the 
excitation of the methylidyne cation CH + and on the intensities of its rotational transitions recently detected in emission in dense 
Q " photodissociation regions (PDRs) and in planetary nebulae. 

Methods. We have developed a nonlocal thermodynamic equilibrium (non-LTE) excitation model that includes the entire energy 
structure of CH + , i.e. taking into account the pumping of its vibrational and bound and unbound electronic states by near-infrared and 
£C) ■ optical photons. The model includes the theoretical cross-sections of nonreactive collisions with H, H2, He, and e~, and a Boltzmann 

£\J 1 distribution is used to describe the probability of populating the excited levels of CH + during its chemical formation by hydrogenation 

of C + . To confirm our results we also performed an extensive analytical study, which we use to predict the main excitation process of 
several diatomic molecules, namely HF, HC1, SiO, CS, and CO. 

■ Results. At densities «h = 10 4 cm -3 , the excitation of the rotational levels of CH + is dominated by the radiative pumping of its 

O electronic, vibrational, and rotational states if the intensities of the radiation field at ~ 0.4, ~ 4, and ~ 300 fim are stronger than 10 5 , 
10 8 , and 10 4 times those of the local interstellar radiation field (ISRF). Below these values, the chemical pumping is the dominant 
source of excitation of the J > 1 levels, even at high kinetic temperatures (~ 1000 K). The far-infrared emission lines of CH + 
observed in the Orion Bar and the NGC 7027 PDRs are consistent with the predictions of our excitation model assuming an incident 
far-ultraviolet (FUV) radiation field of ~ 3 x 10 4 (in Draine's unit) and densities of ~ 5 x 10 4 and ~ 2 x 10 5 cirT 3 . In the case of NGC 
7027, the estimate of the density is 10 to 100 times lower than those deduced by traditional excitation codes. Applying our model 
to other X [ l, + ground state diatomic molecules, we find that HF, and SiO and HC1 are the species the most sensitive to the radiative 
pumping of their vibrational and bound electronic states. In both cases, the minimal near-infrared and optical/UV radiation field 
intensities required to modify their rotational level populations are ~ 10 3 times those of the local ISRF at densities «h = 10 4 cm 3 . 
All these results point towards interstellar media with densities lower than previously established and cast doubts on the dumpiness 
of well-studied molecular clouds. 

Key words. Line formation - molecular processes - radiative transfer - ISM: molecules 

OO 

<N 

l^-v 1. Introduction (PAHs) to pump the molecules to an excited vibrational state 

in bright PDRs (IScoville et al.lll980t iGiardet alJfl997h . While 

■ Our knowledge of the physical conditions of interstellar and cir- this mechanism is neglig ible for certain trans itions (e.g. HF 
-, 1 cumstellar molecular clouds has grown through the confronta- (1-0) and HC1 (1 - 0), Agu ndez et all 1201 lb . the near- and 

\ tion of the observed intensities and shapes of atomic and molec- mid-infrared pumping of the vibrational states has been found 

• • ■ ular lines with those predicted by radiative transfer models, to dominate the rotational excitation of several other species 

. ^ I P uzz l e consists of simultaneously reproducing the emis- (e.g. H 2 0, HNC, NH 3 ) in the circumstellar envelopes of AGB 

k>( ! sion profiles of various species observed in a single astro- stars such as IRC 1021 6 (|Gonzalez- Alfonso &Cernicharoll 19991: 

^ ■ nomical source with the minimal number of physical compo- iGonzalez-Alfonso et al.ll2007t lAgundez et al.ll2008l) . Similarly. 

& \ nents with different temperatures, densities, and velocity fields, the pumping of the electronic states by absorption of UV and 

However, the reliability of the results strongly depends on the optical photons followed by radiative decay has been proposed 

approximations of the radiative transfer models, and in panic- to play a major role in the rotati onal excitation of diatomic 

ular, on the processes expected to drive the di stribution of a molecules dCernicharo et al.l fl997l) . For instance, this process 

mole cule among its excited leve ls. For instance, iMorrisI d 19751) has been studied to account for the near-infrared emission lines 

and ICarroll & Goldsmith! dl98ll) have shown that considering of H? detected in the interstellar and circumstell ar media (e.g. 

only the excitation by nonelastic collisions and the radiative tran- iBlack & van Dishoecl [T987l | Haba rt et al. 1 12005b . the rovibra- 

sitions between the rotational levels of the ground vibrational tional emission lines of CO detected in mole c ular clouds and 

state may lead to an overestimation (by several orders of mag- circum stellar disks (e.g. iKrotkov et al.l Il980t fTroutman et al.l 

nitude) of the density of the main collisional partner or of the 1201 ll). the rotational fine structur e absorption lines of interstellar 



column density of the molecule. C2 dvan Dishoeck & Blackl ll982). the 405 nm band of interstel- 

_ , . r. lar C3 observed in absorption towards HD210121 (Roueffet al. 

Previous theoretical investigations of alternative excitation ~ nn ~? , . ,. C ^ T1 + j I TT 7 — IT - ^ 

. . , . , j. ° , , ■ _c 2002), and the rovibronic emission lines of CH detected in the 



mechanisms have mainly focused on the capacity of near- tz — „ , „ , _ . , , , 1=—. — — - — ii . nn> i 

. , , . . , . , blue spectrum of the Red Triangle nebula Balm & Jura 1992). 

infrared photons emitted by polycyclic hydrocarbon molecules 1 
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Lastly, state-to-state chemistry, i.e. chemical reactions between 
the specific quantum states of the reagents and the products, 
may be involved in the excitation of a molecular species if 
the timescale of the reaction is comparable to that of nonre- 
active and nonelastic collisions (BlacM fl998l) . Indeed, molecu- 
lar excitation via the chemical formation itself has already been 
propos ed to populate the highest rotational levels of interstel- 
lar C3 ([Rp ueff et aL|| 2002l) . and has recently been invoked by 
iGoicoechea et alJ d201 lb as an alternative solution to explain the 
extended and intense rotational emission of OH detected in the 
Orion Bar PDR with the Herschel/PACS instrument. 

Because all these processes could alter our concept of many 
astrophysical environments (from interstellar and circumstel- 
lar PDRs to protoplanetary disks), it is essential to methodi- 
cally explore their respective role on the rotational emission 
lines of each molecular species. While other diatomic molecules, 
such as OH and SiO, will be treated in forthcoming papers, 
we focus here on the methylidyne cation . Since the discov- 
ery of CH + in the interstellar medium by iDouglas & Herzbergl 
( 1 194 lb . the energy structure, spectroscopic properties, nonreac- 
tive collisional processes, and the chemistry of this molecular 
ion have been t he subject of both experimental and theoreti- 
cal studies (e.g.lHakalla etal.ll2006l: lMulledl2010l:lTurpin etaII 
l2010tlGerrich et alJl987tlHierl et alJl997f and reference fherin). 
We have collated this information here into a concise paper with 
the aim to present a complete model of CH + rotational excitation 
and propose an interpretation of its far-infrared emission lines re- 
cently observed with ISO and the Herschel Space Observatory. 



2. Pumping mechanisms of CH 4 

2.1. Energy structure 



Table 1. Equilibrium molecul ar constants (in cm ) of the two 
first electronic states of CH + (Hakall a et al 2006). Numbers in 
parenthesis are power of 10. 



prior to Pears on & Drouinl (120061) . who reported the first lab- 
oratory detection of the CH + 7=1-0 line, all CH + spec- 
troscopic data were deduced from the measurements of the 
rovibronic bands X'£ + - A'n, 'A -' IT and 3 S - 3 II (e.g . 
ICarrel [l968t ICarrington & Ramsavl fl98l lHakallaet all [2006) 
and f rom the pure rotational spectrum observed towards NGC 
7027 dCernicharoetal.1119971) . 

Referring to Fig. [TJ which displays the energy structure of 
the Douglas-Herzberg band system (A'n - X'Z + ), the first X'2 + 
and A'n electronic states of CH + are split into a set of rovi- 
brational levels. Their respective energies E are des cribed by the 
direct matrix elements of the effective Hamiltonian dBrown et alJ 
1979). Taking into account the spin-orbit coupling that results in 
splitting each rotatio nal level in closed- spaced pairs of opposite 
parity (e and / levels. iBrown et al 1fl975l) . the energies of all lev- 
els can be computed as 

E(A, v, 7) = 

T e (A) + co e {A){v + i) - u e x e (A)(v + ±) 2 + a> e y e (A)(v + ±) 3 + 
[B e (A) - a e (A)(v + i) + y e (A)(v + \) 2 + s e (A)(v + ±) 3 ] x 
(7(7 +1)- A 2 ) - 

[D e (A) +p e {A){v +\) + 5 e (A)(v + i) 2 ] (7(7 + 1) - A 2 ) 2 + 
[H e (A) - a He (A)(v + i)] (7(7 + 1) - A 2 ) 3 ± 
i[^(A) + Qr, e (A)(u+i)] (7(7+1)) + 
i ?Di ,(A)(7(7+l)) 2 , 
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Direct spectroscopic measurements of the pure rotational 
spectrum of CH + are both chemically and technically difficult 
to achieve because (a) CH + is an extremely reactive species 
and thus hard to maintain in detectable amounts, and (b) its ro- 
tational transitions lie at submillimeter and far-infrared wave- 
lengths where the sensitivit y of laboratory spectrom eters has 
been limited until recently dPearson & Drouinl I2006I) . Hence, 



Figure 1. Energies (in K) of the ro vibrational levels of the first 
X'E + and A'n electronic states of CH + . u-, 7-, and /"-values 
are given on the left, below, and on the right of each level. All 
allowed radiative transitions are displayed, including pure rota- 
tional (in blue), rovibrational (in purple), and electronic (in red) 
transitions. 
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(1) 2.3. Nonreactive nonelastic collisional processes 



where T e (A), o> e (A), etc. are the equilibrium molecular constants 
of the electronic state A, and where + and - refer to the e and 
/ sublevels. Each level is therefore defined by its total orbital 
angular momentum A, vibrational momentum u, total angular 
momentum 7, an d parity P. With the e quilibrium molecular con- 
stants derived by Hakal kTet al.l (120061) and given in Table Q] the 
electronic, pure vibrational, and pure rotational level structures 
of CH + are significantly separated in energy from one another as 
long as v < 10 and 7 < 13. 

2.2. Nondissociative radiative processes 

In astrophysical environments, the level populations of a 
molecule are controlled by the spontenous and induced radia- 
tive processes and the collisional excitation. As shown in Fig.Q] 
we consider here all possible radiative transitions of CH + that 
meet the selection rules AA = 0, ±1, AS = 0, AJ = (if AA ^ 
0), ±1, and P'P" = -1. The pure rotational transitions lie in 
the mid- and far-infrared part of the electromagnetic spectrum 
from 10 to 1000 //m, the rovibrational lines occupy the near- and 
mid-infrared domain between 1 and 10 /mi, and the rovibronic 
transitions involve ultraviolet and optical photons between 0.3 

and 0.8 /mi. 

With the conventions adopted bv lLarssonl 0983), the absorp- 
tion line oscillator strength / and the Einstein spontaneous emis- 
sion coefficient A of a transition A',u',J' — > A",v",J" of a 
diatomic molecule are given by 
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where m e and e are the mass and the charge of the electron, h 
is the Planck constant, c is the speed of light, and v, M, and S 
are the frequenc y, the transition dipole moment, and the Honl- 
London factors (Herzberg 1950). Under the Born-Oppenheimer 
and the r-centroid approximations, the electronic and nuclear 
motions are separated and M can be expressed as 

M lv'A"v" = qA-v>A»v»R 2 e (rv>v"), (4) 

where r is the r-centroid of the v' - v" vibrational band, q is the 
Franck-Condon factor, and R e is the electronic transition mo- 
ment. 

Colle cting the vibrationa l band oscillator s trengths com- 
puted by lElander et al.l (1 19771) and measured by Wesela k~et alJ 
(2009), a nd using the Franck- Condon factors and r-centroids ob- 
tained bv lHakalla etail ([2006), we found that the electronic tran- 
sition moment of the A l Tl - system has a clearly defined 
functional form (with a correlation coefficient of 0.97) 



R e (r V 'v") = 2.01 x exp(-1.66 r„^«) ea . 



(5) 



Conversely, since we were unable to extract a simple relation 
betwee n the transition dipole momen t and the internuclear sep- 
aration dChackerian & Tipping|[T982l) for the - X l T. + sys- 
tem, we extrapolated Mq ,/o„" from the values com puted at lower 
v' bv lCheng etail (120071) and lFollmeg et alJ dl987l) assuming an 
exponential dependence on v' - v" . Table [2] contains the corre- 
sponding Einstein spontaneous emission coefficients computed 
for the /', /" = 1,0 transitions of the first vibrational bands of 
the A'n - and - X l Z + systems. 



Since we aim to study the CH + excitation mechanisms over a 
wide domain of parameters (i.e. describing different astronom- 
ical sources), and, for example, at kinetic temperature 7k close 
to the equivalent temperatures of the rovibrational transitions 
(> 1000 K), we include all nonreactive collisional excitation and 
de-excitation processes, even between two different rovibronic 
levels, and considering H, H2, He, and e~ as collision partners. 
Unfortunately, there are no experimental data available for the 
CH + collisional rates, and the theoretical data are usually lim- 
ited to the pure rotational levels. 

For collisions with electrons, w e adopted the rate functions 
of T K obtained bv lLim et all d 19991) for v = and 7 up to 7. For 
collisions with neutrals, we us e d the CH + -H e col lisional rates 
compu ted bv lHamm ami et al. (2008, 2009]) and iTurpin et al.l 
(2010) for u — and 7 up to 10, and then applied the rigid rota- 
tor approximatiorQ to take into account collisions with H and H2- 
For 30 0K < T K < 2000K, we used the data of lHammami et al.l 
(2009). Below this temp erature range, we pr eferentially adopted 
the values computed bylTurpin et al.l ([2010) rather than those of 
lHammami et al.l d2008l) because the former authors have found 
that including of the long-range charge-induced dipole potential 
of CH + has an impact on the collisional rates at low tempera- 
ture. Finally, for 7k higher than 2000 K, we extrapolated the 
de-excitation rates by collisions with neutrals from the values 
computed at lower temperature assuming a power law depen- 
dence on Tk- 

These nonreactive collisional rates were used for all tran- 
sitions with A' = A" and v' = v". In turn, the v' = v" + 
1 collisional rates were deduced from the theoretical cross- 
se ctions for inelastic collision s between CH + and He computed 
bv lStoecklin & Voroninl d2008l) for v' , v" — 1,0 and J', J" up to 
0, 5. For all other transitions occuring between two different vi- 
brational or electronic levels, the downward rates were inferred 
from those of the pure rotational transitions assuming a power- 
law dependence on (E' - E"), 



K A'v'J' 



A"v"J" 



(27' + 1)(27" + 1) X a 



E' - E' 



(6) 



where E' and E" are the energies of the upper and lower lev- 
els and a and b are the best-fit coefficients of all available data 
at the temperature 7k- As a justification for this choice, the 
obtained scaling relations are found to describe the pure rota- 
tional transition rates within a factor of 3 and with a corre- 
lation coefficient ranging from 0.8 to 0.9. We note, however, 
that they systematically underestimate the vibrat ional transition 
rates e stimated from the cross sections of IStoecklin & Voroninl 
(2008) by a factor varying from 70 to 5 for 10 < T K < 70 
K, because the vibrational and rotational quenching cross sec- 
tions have the same order of magnitud e in the low collision en- 
ergy limit dStoecklin & Voroninl l2008). and by a factor < 5 for 
7k > 70 K. Nonetheless, since the discrepancies diminish when 
7k increases and because our subsequent analysis and predic- 
tions will mainly come from the high kinetic temperature do- 
main, we consider our fitting approach as valid. 



1 Within the rigid rotator approximation, the collisional rates are pro- 
portional to yu~ 1/2 , where yu is the relative mass of the system collider- 
CH + . 
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Table 2. E instein spontaneous emission coefficients Aav/'avv ( m s ) computed for the J', J" — 1, tra nsitions. The data are 
taken from lHakalla et alJ (l2006h for the CH + A : n - band system (top) and from lCheng et alJ (120071) and lFollmeg et alJ (Il987h 
for the CH + - band system (bottom). Numbers in parenthesis are power of 10. 
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" Extrapolated values assuming Mo,/oo = 0.65 x exp[-4.61 v'] eao, and Mvoi = 0.50 x exp[-3.38 (i/ - 1)] eao. 



2.4. Reactive collisional processes and photodissociation 

Finally, since the chemical timescales of CH + are com- 
parable with those of the nonreactive nonelastic collisions 
dBruderer et al.ll20Tol) . we took into account the excitation and 
de-excitation of CH + during its chemical formation and destruc- 
tion. 

During the past decades, the advances of the cross molecu- 
lar beam experiments and of the theoretical studie s of chemical 
reaction dynamics (see the reviews of the field by ICasavecchial 
200(i lLevinell2005l) have led to measurements and calculations 
of state-specific reaction rates for several neutral -neutral reac- 
tions (e.g. Q( 3 P) + H 2 (f = l)^OH(u= D + Hbv lBalakrishnanl 
120041) . For a few reactions these studies have even provided de- 
tailed diagrams of the most probable state of the products de- 
pending on the reag ent excitation state and the collision energy 
(Cas avecchial [2000 and references therein). Similarly, the re- 
cent improvements in t he flowing afterglow app aratus and the 
ion trapping techniques dGerlich & Bor odi 2009) have led to the 
first experimental deductions of state-specific rate coefficients of 
C + + H2 <^ CH + + H: depending on the internal energ y of H2 
(up to v = 1) for the forward reaction (Hier l et aUl 19971) and on 
the internal energy of CH + (up t o J = 3) for the rev erse process 
dPlasil et al.ll2011l) . In particular. IPlasil et alJ d201 lb showed that 
the destruction rate of CH + (7 = 0) by reactive collisions with H 
is 2 to 30 times lower than that of rotationally excited CH + . 

To include the CH + st ate-to- s tate c hemi stry, we followed 
the ap proach described by iBlackl (1 1998b and Ivan der Tak et alJ 
J2007h . 

• We considered that CH + (A, v, J) is mainly destroyed by re- 
actions with H, H2, and e~ or by photodissociation. Because 
of lack of information, we also assumed that the three last 
reactions are independent of A, v, and J by adopting generic 
reaction rates for all levels: IcR = 1.2 x 10~ 9 cm 3 s , kP- — 

ill e 

1.5 x 10" 7 (r D /300Kr ' 42 cm 3 " s _1 , and*£ = J %(r y (v)I v dv 
s~l for reactions wit h H2 and e~ and for the photodisso- 
ciation respectively dMcEwan et al.l 119991: iMitchelll Il990t 
Ivan Dishoeck et aLll2006l) . where 7Yi is the temperature of 
the reaction and <x r (v) are the CH + photodissociation cross- 



sections computed by iKirbvetaD (Il980j). Convers ely, we 
adopted the state-specific rates of IPlasil et alJ (1201 ll) for the 
destruction of CH + by collisions with H atoms. 
• During the chemical formation of CH + we assumed that the 
rotational level populations follow a Boltzmann distribution 
at the temperature of the reaction Tp: 

k iv,j = kF + 1) X exp(-£ w /r F ) / Q(T F ), (7) 

where k F is the state-averaged formation rate and Q(Tp) is 
the partition function of CH + . 

The formation and destruction temperatures correspond to the 
collision energies of the reagents and thus depend on their in- 
dividual dynamics and on the kinetic temperature of the gas. 
However, since our subsequent analysis is based on the results 
of static models at chemical equilibrium, we assumed in the fol- 
lowing Tp, = Tp = Tk, unless specified otherwise. Because the 
dynamical properties of the gas are neglected, this hypothesis is 
a lower limit on the actual values of I'd and Tp. 

2.5. Critical parameter analysis 

In practice, all the above processes are not necessarily signifi- 
cant over the entire domain of parameters, namely the density 
of the gas «h (in cm -3 ), the kinetic temperature 7k (in K), 
and the intensity of the radiation field at all wavelengths I v (in 
erg crrT 2 s~' Hz 1 sr~'). Moreover, if they are all included in 
the radiative transfer models, they can lead to computational 
errors when the populations of the high-energy levels are too 
small. It is therefore useful to make a preliminary exploration 
of the parameter domain and reduce, if necessary, the number 
of levels and excitation processes considered depending on the 
physical conditions of the gas. To do so, we used the analytical 
treatment of the main de-excitation mechanisms that we present 
in Appendix lAl In addition to the critical density «h cri.i, we de- 
rived for each level i of a molecule four dimensionless critical 
parameters: X CI [j, Y cr [j, Z CI [j, and S cr ij which depend on the re- 
active and nonreactive collisional rates, the abundances relative 
to H of the collision partners, the Einstein spontaneous emission 
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coefficients, and the intensities of the radiation field at the corre- 
sponding wavelengths. 

Given the energies and the lifetimes of the A > 1 or v > 1 
levels, collisional de-excitation of vibrational and electronic lev- 
els are found to be always negligible compared to radiative de- 
excitation. More importantly, our approach shows that the pa- 
rameter space is divided for each rotational level into six regions 
where the de-excitation occurs mainly by (1) reactive collisions, 
(2) nonreactive collisions, (3) pure rotational, (4) rovibrational, 
(5) rovibronic radiative transitions, and (6) photodissociation. 
For example, if «(H2)/«h = [H2] = 0.5 and n(e")/ n H = [<?"] 
10 ', the limits of the first five regions for CH + in the level 
K,vJ = 0,0, 1 are T K ~ 800 K, n H ~ 10 7 cm 4 , 7„(4 /im) ~ 
2 x 10" 8 erg cm" 2 s" 1 Hz" 1 sr" 1 , and 7 y (0.4 /an) ~ 7 x 10" 12 
erg cm" 2 s" 1 Hz" 1 sr" 1 . 

Since the excitation processes are not included, the critical 
parameter approach is insufficient to calculate the distribution 
of a molecule among its levels. However, it proves useful to 
point out the physical conditions for which several excitation 
and de-excitation mechanisms have to be taken into account and 
the level populations must be computed with non-LTE radiative 
transfer models. 
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Figure 2. Interstellar radiation field intensity I v at high Galactic 
latitude from Tableland Eqs.[8]-[T4]as a function of the wave- 
length. The wavelengths of the 1725 allowed pure rotational, 
rovibrational, and rovibronic radiative transitions of CH + are 
also displayed in blue, purple, and red, taking into account the 
first two, five, and thirteen electronic, vibrational, and rotational 
levels ofCH + . 



manner as 



3. Radiative transfer modelling 

3.1. Method: the MADEX code 



Table 3. Reference parameters used in Eqs.|8]-[T4]to compute 
the interstellar radiation field at high Galactic latitude. Numbers 
in parenthesis are power of 10. 



frequency domain 


index 


T(K) 


T 


v (GHz) 


P 


microwave 


CMB 


2.7 








far-infrared 


fir 


18 


1.7 (-05) 


1.1(3) 


2.0 


mid-infrared 1 


mi[ 


50 


2.5 (-08) 


2.9 (3) 


1.0 


mid-infrared 2 


mii 


260 


3.0 (-10) 


1.5 (4) 


1.5 


near-infrared 


nir 


3000 


6.0 (-13) 


1.8(5) 


1.0 



To study the individual impact of all pumping processes 
and compute the steady-state level populations of CH + in sev- 
eral astrophysical environments, we improved and updated the 
MADEX ( MADrid molecula r spectroscopy EXcitation) excita- 
tion code dCernicharol |20 1 2|) . MADEX is a radiative transfer 
model bas ed on the multi-shell high velocity gradients (LVG) 
formalism (Goldrei ch & Kwanll 19741) . The physical model con- 
sists of a spherical (or plan parallel) cloud of gas and dust defined 
by the radial (or perpendicular) profiles of the density the H2 
He, and e~ abundances relative to H, [He], [H2], and [e ]; the 
CH + column density A^(CH + ); the kinetic temperature T^; the 
expansion velocity of the cloud along the line of sight f exp ; and 
the dust abundance and temperature and Tj. 

Gas and dust are pervaded by an internal radiation field emit- 
ted by a central star with a radius and an effective temperature 
r*, and an external radiation field defined by its specific inten- 
sity I v . To cover a broad range of the electromagnetic spectrum, 
the external radiation field is modelled in a simple and versatile 



/v = 



2/z v 3 


1 


c 2 


exp(/!v/fcr CM B) - 1 


2/1 v 3 




c 2 


exp(/iv/fc7fi r ) - 1 L 


2/1 v 3 




c 2 


exp(hv/kT mil ) - 1 


2h v 3 


Xmi 2 


c 2 


exp(hv/kT mil ) - 1 


2/zv 3 


Xnii [ 



(8) 



r —— [l-expC-Tfcfv/vfl^)] (9) 

1 fir) - 1 L J 



— [l-exp(-T mi] (v/v rai /».)] (10) 



- [l-eXpC-T^Cv/Vmi,^)] (11) 



- [l - expC-r^Cv/v^)] (12) 



c 2 exp(hv/kT nh ) ■ 
+ Xo P t 9.4 10 6 v" 1 - 7 if 3.7 x 10 14 < v< 1.5 x 10 15 GHz (13) 



Xuv 

4tt 



6.36 10~ 17 v 2 



1.02 10" 21 " 3 



v 4.08 10~ 27 v 4 

— H : 



. (14) 



where Eqs. [8] - [14] correspond to the cosmic microwave back- 
ground emission, the emissions of the three dust components, 
the large amorphous carbon and silicates, the small amorphous 
carbon grains, and the PAHfH as well as the emission of cold, 
A-type, and OB-type stars. As a reference, the temperatures and 
opacity laws of these functions were chosen to match the in- 
terstell ar radiation f i eld at high Galac tic latitud e in the ultra- 
violet iDrainel 1 19781; [Henry et all 1 19801). optical jvan Dishoeckl 
1 1 9881: iKoppl 19961). near^nfrared (Moskalenko et al.ll2006l) . mid- 
infrared ( Are ndt et al.l Il998t ICom piegne et al. 2011), and far- 
infrared dLagache et al.|[l999t iFinkbeiner et al.l ll999) domains. 
Each component can also be adjusted using the scaling factors 
ATfir. #mi, . Afmi 2 , *nir, ATopt, and Xuv The reference parameters are 
given in Table [3] and the corresponding specific intensity of the 
local interstellar radiation field is displayed in Fig. [2] along with 
the wavelengths of the pure rotational, rovibrational, and rovi- 
bronic transitions of CH + . 



2 for the sake of simplicity, a smooth e mission funct i on is 
adopted, neglecting t he PAHs band emitters ([Draine ~&T1 l200l 
ICompiegne et al 1201 lh . 
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Lastly, the absorption of the internal and external radiation 
fields at near-infrared, optical, and ultraviolet wavelengths by 
dust particles across the cloud is taken into a ccount in MADEX 
assuming the extinction curve derived by Fitzrjatrick ( 1999) 
and th e relation A v = Nn (cm"~ 2 )/1.8 x 10 21 dPredehl & Sc hmitt 
fl99l between the visual extinction Ay and the total column den- 
sity N H - 

3.2. Analysis of the main excitation pathways 
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Figure 3. CH + steady-state normalized excitation diagrams 
computed with MADEX assuming «h = 10 4 cm" 3 , 7k = 50 
K, Xo? t = 10 5 , and N(CU + ) = 10 13 cm 2 . The number of vi- 
brational levels taken into account within the X l "L + and A l Tl 
electronic states varies between one (blue squares) and five (red 
pentagons). We also display for comparison the predictions of 
MADEX if the radiative pumping processes are switched off 
(black crosses). 



Our first goal was to find the optical, near-infrared, and far- 
infrared radiation fields required to perturb the steady-state rota- 
tional level populations from what they would be in the absence 
of radiative pumping. Therefore, while more advanced models 
are discussed in the next section, we consider here a spherical 
cloud of gas only (no extinction from dust) without a central star 
and with constant density, temperature, and chemical composi- 
tion. Moreover, to focus on the radiative pumping of bound rovi- 
bronic states, we intentionally neglected in a first step the pho- 
todissociation process. For [H2] = 0.5, [He] = 0.1, [e~] = 10" 4 , 
vV(CH + ) = 10 13 cm" 2 , and v exp = 1 km s" 1 , we computed 75600 
models (for a total computing time of ~ 19 hours): 

• three densities: «h = 10 3 x 10 !_I for i = 1 to 3; 

• 60 kinetic temperatures: 7"k = 30 x 1.06' -1 for i = 1 to 60; 

• 140 far-infrared, near-infrared, and optical radiation fields: 
Xat = 10 3 x 1.09'" 1 for i = 1 to 140,*,* = 1, and* opt = 1, 

= 1, Afrit = 10 6 x 1.09'- 1 for i = 1 to 140, and^- opt = 1, 
and^ fir = l,Xnh- = 1, and^ p t = 10 2 x 1 .1 l f 1 for i = 1 to 
140, 

i.e. exploring the parameter space below the limits set by the 
critical parameters analysis (see Sect. l2.5l and AppendixlAl. The 
dominant excitation and de-excitation mechanisms are discussed 
in Appendix FBI and are shown in Fig. IB. II As expected for this 
parameter range, the de-excitation of the rotational levels mainly 
occurs via the pure rotational radiative transitions for all models. 
Conversely, Fig. lB.ll reveals five different excitation regimes. 



1 . In the limit of low kinetic temperature, the reactive collisions 
are responsible for the excitation of the / > 1 rotational lev- 
els, implying «ooy K nnexp(-Eooj /T^) (see Eq. IB. Il l at the 
steady-state equilibrium. 

2. In the limit of high kinetic temperature, the J > 1 rota- 
tional levels are excited through chemical pumping of a high- 
energy rotational state (e.g. 7=10) followed by the radiative 
decay of all intermediate J states. Interestingly, the nonreac- 
tive collisions contribute to less than 20 % of the rotational 
excitation even for high values of 

3. In the limit of large far-infrared radiation field, the rota- 
tional excitation occurs through the successive excitations of 
all intermediate rotational levels by direct absorption of far- 
infrared photons. At the steady-state equilibrium, we derive 
"00/ ^A'fir ( see Eq.|B3]l. 

4. In the limit of large near-infrared radiation field, the pure ro- 
tational levels are excited by near-infrared pumping of the 
first vibrational level of the X l T. + state. Since the transition 
dipole moments Mo„oo (for v > 1) are at least two orders 
of magnitude smaller than M0100 (see Tabled, the pumping 
of the v > 1 levels has a negligible contribution. Moreover, 
because the transitions are restricted to AJ = +1, the den- 
sity of CH + in the level J results from the successive exci- 
tations of all intermediate levels. Consequently, the high-/ 
levels are more difficult to populate and are more sensitive 
to the intensity of the near-infrared radiation field: at steady- 
state nooj oc^ +1)/2 and n mj ^ ( se e Eqs.lB.4l&lB3Tl for 
odd and even values of J. 

5. Finally, in the limit of large optical radiation field, the pure 
rotational levels are excited via the optical pumping of the 
first rotational levels of the A 'II vibronic states followed by 
the radiative decay through the rovibrational levels of the 
X l I. + state. As shown in Fig. IB. II the selection rules imply 
that all vibrational levels u < J - 2 are required to accurately 
compute the population of the pure rotational level J. This 
result is illustrated in Fig. [3] which displays the normalized 
excitation diagrams of CH + (up to J — 6) computed with 
MADEX for n H = 10 4 cirT 3 , T K = 50 K, X o P t = 10 5 and 
taking into account an increasing number of vibrational lev- 
els (from 1 to 5). Lastly, since the excitations of all pure ro- 
tational levels ultimately originate from the same transition 
A,v,J- 1 , v, 1 <— 0, 0, 0, their populations are equally sen- 
sitive to the intensity of the optical radiation field: «oo/ rc Xopt 
(see Eq. IB.6I > for all J > 1 at the steady-state. 

3.3. Radiation field required to activate the radiative pumping 

With our grid of models we derived the radiation field intensities 
above which the far-infrared, near-infrared, and optical pump- 
ings play a dominant role in populating the excited levels of 
CH + . These limits are displayed in Fig. [4] for the J < 5 rota- 
tional levels as functions of the kinetic temperature and the den- 
sity. The main excitation processes are found to exponentially 
depend on Tk, as was already predicted in the previous section. 
In addition, the superimposition of the curves computed at dif- 
ferent densities confirms their dependence on «h and / derived 
analytically in Appendix iBl 

Fig. |4] shows that the interaction with the far-infrared radi- 
ation competes with the collisional processes in the rotational 
excitation of CH + for Xfr > 10 4 x (n H /10 4 ) l/J . Similarly, the 
pumping of the vibrational states could influence the popula- 
tions of the low-J rotational levels, providing that the gas is per- 
vaded by a near-infrared radiation field at least 10 7 x («h/10 4 ) 2/ ' / 
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Figure 4. Far-infrared (top), near-infrared (middle), and optical 
(bottom) radiation fields (in units of the local ISRF) at which 
the pure rotational levels of CH + are equally populated by colli- 
sions and radiative pumping, depending on the kinetic tempera- 
ture and the density of the gas. Below these curves, excitations 
are primarily driven by reactive collisions and chemical pump- 
ing. Above these curves, excitations are primarily driven by di- 
rect radiative pumping (top) and by the radiative pumping of the 
vibrational levels of the X'2 + state (middle) and of the A'll elec- 
tronic state (bottom). 



and 10 s x (n H /10 4 ) 2/y times larger than the local ISRF for 
T K < 100 K and T K > 100 K. Lastly, an optical radiation field 
of ~ 10 5 x («h/10 4 ) the local ISRF is required to activate the 
pumping of bound electronic states for 7k > 100 K, while this 
limit drops at lower kinetic temperatures. 

3.4. Impact of the photodissociation 

As emphasized in AppendixlAl the impact of the photodissocia- 
tion on the rotational excitation of CH + strongly depends on the 
spectrum of the UV radiation field. For instance, while the pho- 
todissociation is negligible compared to bound-bound electronic 
transitions in a radiation field emitted by a black body with an 
effective temperature < 10 4 K, the two processes have compa- 
rable rates in the standard interstellar radiation field. Therefore 
the photodissociation could play an important role in the rota- 



tional excitation of CH + in any environment in which the A'll 
- X l Y, + pumping competes with the collisional processes (Fig. |4] 
bottom). 

Hereafter, we therefore included the photodissociation pro- 
cess into MADEX and applied a model to several astronomical 
sources with different values of Tr, Xfa I '«h, Xmi Inn, and^ opt /nH- 
The following questions were addressed: how do the predictions 
of MADEX compare with the available observations? How do 
all these excitation processes modify our understanding of the 
physical conditions of the interstellar and circumstellar matter? 

4. Application to astronomical sources 

4.1. Chemical formation of CH + 

Because an efficient production pathway is required to balance 
its rapid destruction by hydrogenation and dissociative recombi- 
nation, CH + is believed to be formed either via 



H? 



CH + + H + 



(15) 



in environments submitted to a strong X-r ay radiation field (e.g. 
the central molecular zone, lLanger|[l978l) . or by the highly en- 
dothermic reaction 



C + + H 2 ^ CH + + H (AE/k = 4640/Q 



(16) 



in other Galactic environments. So far, two main scenarios have 
been invoked to overcome the endothermicity of reaction[16] (1) 
by the release of supra-thermal energy, induced for example by 
magnetohydrodynamic shocks, in regio ns with large reservoirs 
of kinetic and magnetic en ergies (e.g. lDraind[l986t IXie et al] 
11995 ; Falgarone et al.ll 1995h; and (2) by t he internal energy of 
Ht dHierl et aljll997HAgundez et all ZOlO) in regions where the 
vibrational levels of H2 are excited by FUV fluorescence. 

For each astronomical source considered below, we there- 
fore applied an appropriate chemical model to describe the abun- 
dance distributions of CH + and of its main collisional partners. 
To investigate the origine of the rotational excitation of CH + 
we then systematically ran MADEX in three different config- 
urations: O considering only the excitation by nonreactive colli- 
sions, © considering all collisional excitation processes (i.e. in- 
cluding the chemical pumping), and © considering all excitation 
processes (i.e. including the radiative pumping and the photodis- 
sociation). The parameters and results of all models are summa- 
rized in Table [4] where we compare the continuum-substracted 
line fluxes (or continuum-substracted line intensities) of CH + 
computed with MADEX with the available observations. 

4.2. Diffuse interstellar medium 

We first focus on the cold diffuse interstellar matter character- 
ized by a low density and by low^^/we and^ opt /nH ratios. The 
observed mean column density of CH + , log in (A^( CH + )/A^H) ~ 




0.33 
Gredell 



(from a compilation of d a ta: Crane et al 



19971: IWeselak et al.l l2008t [Falgarone et al 
. iGodardet al.l l2012h has been proposed to trace low- 
veloc ity magnetohydro dynamic shocks dPineau des Forets et al .1 
1986 ), Alfven waves ( Federman et al.l Il996h . tu rbulent mix- 
ing (jLesaffre et al.1l2007 ), or turbulent dissipation dJoulain et al] 
1998). According to t he TD R (turbulent dissipation region) 
model of IGodardet "all (2009), CH + is formed in regions where 
turbulent dissipation is triggered by the decoupling of the ion- 
ized and neutral flows at small spatial scales (~ 100 AU). In 
these regions, the kinetic temperature increases to ~ 1000 K, 
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Table 4. Comparison between the observed intensities J Idv or fluxes JfldvdQ. of the rotational lines of CH + and those predicted 
by MADEX for the selected astronomical sources. In all cases, MADEX was run in the three configurations O, ©, and © (see main 
text), and the far-infrared continuum emission was removed in the computation of the line intensities. 



parameter 




cold diffuse medium 






Orion Bar 




NGC 7027 


iH(crrr J ) 






50 






5xl0 4 




2x 10 5 


Xopt 






1-50 






3x 10 4 




4x 10 4 


Xnir 






1-50 






3x 10 4 




2 x 10 5 








1 - 10 






1 




1 


T Sl 






18 






55 




150 


rfli 






1.7 x 10~ 5 






5 x 10~ 2 




1 x 10-2 








2.0 






1.75 




0.5 


fexp (km s _1 ) 






5 






10 




20 


fi(sr) 












2.7 x 10~ 8 




1.2 x 10~ 9 


transition 






intensity 






flux 




flux 




(in 10" 


10 erg s -1 cnT 2 sr~') 




(in 10 


-14 -1 -2\ 

erg s cm ) 




(in 10~ 1S W cm" 2 ) 




O 


© 


© obs a 


O 


© 


© obs" 


O 


© © obs r 


7=1-0 


2.7 


4.5 


5.9 < 79 


9.3 


11.8 


11.8 11.1 + 1.8 


0.4 


0.5 0.5 0.47 ± 0.01 


7 = 2-1 


3.7 


7.1 


9.5 < 480 


13.2 


18.1 


18.3 


1.6 


2.0 2.0 1.51 ±0.05 


7 = 3-2 


2.4 


7.0 


9.7 


6.7 


13.8 


13.9 


2.1 


2.5 2.5 2.18 ±0.17 


7 = 4-3 


2.3 


7.7 


1.1 


6.1 


13.4 


13.5 


1.6 


2.1 2.1 2.00 ±0.22 


7 = 5-4 


1.1 


6.6 


9.6 


3.0 


8.3 


8.4 


1.1 


1.7 1.7 2.50 ±0.41 


7 = 6-5 


1.1 


6.3 


9.2 


2.9 


7.5 


7.6 


0.7 


1.3 1.3 2.41 ±0.33 



" dete ction limits computed with the r ms noise level of th e C H + (1-0) and (2- 1) spect ra with the Hersch el/HIFI instrument 1 Goda rdet al.l 
2012; M. Gerin, priv. comm.); * from lNavlor et all j201fj) and lHabart etafl d2010h : c from lCemicharo et all dl997h and l Wesson et al] fcOlCh . 



while the CH + formation temperature is boosted by the ion- 
neutral velocity drift up to ~ 2000 K. Assuming a density 
«h = 50 crrT 3 and a moderate shielding Ay > 0.2 of the in- 
terstellar radiation field (y opt = Xmr - X&r = 1)> most of the hy- 
drogen is in molecular form and the ionization degree is ~ 10 -4 . 

These properties were injected into MADEX to compute the 
beam-averaged continuum-substracted intensities of the CH + ro- 
tational lines (see Table |4]l assuming that the turbulent dissipa- 
tion regions have a surface filling factor of 1 . We find that the 
pumping by near- and far-infrared radiation has no influence on 
the excitation of the rotational levels of CH + even if Xmr = 50 
and Xfc = 10, i.e. adopting the high est energy density of the 
interstellar radiation field derived by Moskalenko et al. (2006) 
across the Galactic disk. In contrast, the optical pumping and 
the photodissociation are found to be non-negligible compared 
to the collisional processe s for y uv / npt > 10, i.e. in the most cen- 
tral regions (R g < 4 Kpc, Moskalenko et al.ll2006h of the Milky 
Way. For instance, for^- uv / opt = 50, the optical pumping and the 
photodissociation contribute up to 30 % of the rotational excita- 
tion of CH + (Col. 4 of Table EJl. In all cases, we finally find that 
the intensities of the two first rotational emission lines predicted 
by MADEX are at least 14 times smaller than the li mit of de- 
tection of the Herschel/HIFI data (God ard et al.ll2012l) . in agree- 
ment with the fact that CH + has only been detected in the diffuse 
interstellar medium in absorption towards strong submillimeter 
continuum sources. 



4.3. Hot and dense PDRs: e.g. the Orion Bar 

We selected the Orion Bar as a prototypical hot and dense PDR 
wit h an high y n nt/wH r atio. Lying at a distance of 414 + 7 
pc dMenten et al] [2007), the Bar consists of a homogeneous 



angular size (") 
5 10 



interclump medium wi th a density «h 



5 x 10 4 crrT 



(I Young Owl et al.l l2000). located approximatel y 111" (~ 0.22 
pc) s outh-east of the Trapezium star cluster (Pellegrin i et al] 
2009). The resulting incident UV and optical radiation fields 




Figure 5. Kinetic temperature and abundances relative to H of 
H2, H, e~, and CH + computed with the Meudon PDR code as 
functions of the extinction (bottom axis), the distance (middle 
axis), and the angular size (top axis) across the Orion Bar starting 
from the ionization front. 



at the ionization fron t are ~ 3 x 10 4 that of the local ISRF 
dMarconi et al 1 fl998h . The near-infrared continuum emission 
of the HII region cor responds to 10 4 that of the local ISRF 
dWalmslev et al] l2000h . Lastly, the rece nt observations a t 70, 
160, 250, 350, and 500 /mi performed bv lArab et ail d2012l) us- 
ing the PACS and SPIRE instruments of the Herschel space tele- 
scope have revealed dust temperature and spectral index gradi- 
ents from 80 K to 40 K and from 1 .2 to 2.2 from the ionization 
front to the Orion molecular cloud, respectively. 

With all these parameters, we computed the chemical profile 
of the Orion Bar with the Meudon PDR code, a one-dimensional 
chemical model in which a static slab of gas of given thickness 
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is illuminated on one side or on b oth sides by a given FUV ra- 
diation field dLe Petit et alJl200 6). Fig. [5] displays the physical 
structure of the Orion Bar predicted by the most recent online 
version of the Meudon PDR codfl which takes into account the 
influence of vibration ally excited H2 in the formation of CH + 
( Agu ndez et alj|20Tol) . Inserting the density profiles of H, H2, 
e~, and CH + into MADE X and taking into accou nt the geometry 
of the bar described bv lHogerheiide et al. ( 1995), we finally de- 
rived the six first rotational line fluxes of CH + (integrated over a 
36" wide solid angle, after having removed the far-infrared con- 
tinuum emission) given in columns 6, 7, and 8 of Table |4] 

The comparison of models O, ©, and © shows that the ra- 
diative pumping and the photodissociation have only a marginal 
effect (< 2 %) on the excitation of the rotational levels of CH + . 
In contrast, even if the high-/ transitions appear to be mostly 
driven by chemical pumping, both reactive and nonreactive col- 
lisions are required to accurately model all the far-infrared emis- 
sion lines of CH + . At last, the flux of the CH + 7=1-0 transition 
derived w ith MADEX i s foun d to agree well with the value ob- 
served bv lNavlor etail (|2010) and Hab art et all (1201 Ol) using the 
Herschel/SPIRE instrument. This result indicates that the emis- 
sion lines of CH + detected in the Orion Bar probably arise from 
the interclump medium rather th an from dense and hot r egions 
(T K ~ 200Kand« H ~ 10 7 cm~ 3 , lGoicoechea et alj201 lb . Since 
the J = 3 - 2 transition of CH + detected with Herschel/PACS 
was found to spatially correlate with the 60-100 fj.m rotational 
lines of OH dGoicoechea et al.l [201 lh and t he high- J rotationa l 
lines of CO detected with Herschel/SPIRE dHabart et alJl2010h . 
this result also suggests that chemical formation and/or radia- 
tive pumping could have a significant impact on the rotational 
excitation of OH and CO in the Orion Bar PDR. If this is the 
case, clumps would no longer be required to account for their 
observed emission line intensities (Goi coechea et alj|201 lh . 



Lastly, we chose NGC 7027 as an archetypal planetary neb- 
ula with an high y „\ r / «h ratio. Located at a distance of 880 + 150 
pc dMassonll 19891) . NGC 7027 is composed of ionized, atomic, 
and molecular C-rich shells surrounding a single central star 
with an effective tem perature T+ ~ 1 98,000 K and a radius 

~ 5.21 x 10 9 cm dLatter et al.ll2000l) . The ionized shell re- 
vealed by a strong near-infrared continuum emission has an el- 
lipsoidal morphology with dimensions 6" x 12", while the bor- 
ders of the PDR traced by the emission line s of H2 have a bi- 
conal structure with dimensions 10" x 13" dLatter et al.l f2000: 
ICox et al1l2002l) . The molecular envelope detected in CO and 
other millimetric molecular transitions extends up to ~ 100" in 
radius dNavarro et al.ll2003l) . 

For the sake of simplicity, we modelled the PDR region 
as a spherical shell with dimensions set to match the molec- 
ular e mission lines o bserved along the minor axis of NGC 
7027 (ICox etal.1 120021) . The ionization front is located at 3.4" 
(~ 0.0145 pc) from the central star, implying an incident UV ra- 
diatio n field of ~ 3 x 10 4 that of the local ISRF dHasegawa et al.l 
2000 ). Using the optica l and near-infrared spectra modelled 
bv IVolk&Kwokl dl997) and the extinction curve found by 
IZhang et all Q2003), we estimate that the incident optical and 
near-infrared radiation fields are ~ 4 x 10 4 and ~ 2 x 10 5 
times stronger than the local ISRF at 0.4 /im and 4 fim. Finally, 
the far- and mid-infrared dust emission are modelled as a sin- 
gle grey body with a temperature T^ T = 150 K, an opacity 
Tfi r = 10 -5 , and a spectral ind ex R^x _ 0-5 to reprod uce the con- 
tinuum emission observed by Dvc k & Simonl (1 19761) and Terzian 
( 1989b . In contrast to t he work of previous au thors dCox et al.l 



2002; lHasegawa etail l2000t lYan et al.l 1 19991) . who proposed 
PDR models including a high-density shell, we considered only 
the low-density PDR component defined by «h = 2 x 10 5 cm -3 . 



4.4. Planetary nebulae: e.g. NGC 7027 
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Figure 6. Kinetic temperature and abundances relative to H of 
H2, H, e~, and CH + computed with the Meudon PDR code as 
functions of the extinction (bottom axis) across the PDR region 
of NGC 7027 starting from the ionization front. The distance and 
the angle from the central star are indicated along the middle and 
top axes. 



version 1.4.4 available at http://pdr.obspm.fr/PDRcode.html 



In the same way as for the Orion Bar PDR, the chemi- 
cal composition of NGC 7027 was computed with the Meudon 
PDR code, i.e. assuming that the radial structure of the source 
is well described with a plan-parallel model. The carbon-rich 
circumstellar envelope was modelled with C, O, N, and S ele- 
mental abu ndances of 1.3 x 1 0~ 3 , 5.5 x 10~ 4 , 1.9 x 10~ 4 , and 
7.9 x 10~ 6 (Middlemass 1990). The resulting chemical and tem- 
perature profiles are shown in Fig.|6]as functions of the distance 
from the central star. The fluxes of the rotational emission lines 
of CH + (integrated over a 9" wide solid angle, after removing 
the far-infrared continuum emission, see Table |4]i were finally 
calculated, assuming that the chemical composition has a spher- 
ical symmetry and taking into account the influence of the limb 
brightening. 

Because of the high kinetic temperature and electronic frac- 
tion, the nonreactive collisions of CH + with e~ pilot the dis- 
tribution of CH + among the low-/ rotational levels, and even 
compete with the chemical pumping in the excitation of the 
high-/ transitions. Our main result is that the predictions of 
MADEX agree excelle ntly with the observations performed by 
ICerni charo et al. (1997) using the ISO/SWS instrument and with 
those of IWesson et al.l yOlO) using the Herschel/SPIRE instru- 
ment. Furthermore, because all collisional excitation processes 
are included in MADEX, the intensities of all observed rota- 
tional lines are explained with a circumstellar medium of den- 
sity «h = 2 x 10 5 , a value 10 to 100 times lower than those 
inferred from traditional excitation model s dCernicharo et al.l 
U997tlHasegawa et alJ l2000: Yan et al.lll999h . 
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4.5. Summary 

We conclude from these applications that the far-infrared pump- 
ing, the near-infrared and optical fluorescences, and the pho- 
todissociation have a marginal impact on the far-infrared emis- 
sion lines of CH + . For all the sources considered, the excitations 
of the high-/ and low-/ lines mainly arise from chemical pump- 
ing and from a combination of chemical pumping and nonreac- 
tive collisional processes. For low-/ lines, the balance of these 
two mechanisms then principally depends on the chemical com- 
position of the gas and the kinetic temperature. For instance, in 
environments with a large electronic fraction, such as the C-rich 
PDRs (e.g. NGC 7027), the nonreactive collisions with electrons 
are found to dominate the rotational excitation of CH + up to the 
/ = 4 — 3 transition. 

Our main result is that the comparison between the observed 
line fluxes and the predictions of MADEX points towards in- 
terstellar and circumstellar media with densities far lower than 
those deduced with traditional models. This study thus proves 
the importance of exploring all possible excitation pathways of 
a molecule before drawing any conclusion on the physical struc- 
ture and geometry of the source. In particular, we question the 
degree of dumpiness of the Orion Bar molecular cloud. 

We therefore recommend to treat any other astronomical 
environments detected in CH + in a similar way. For instance, 
the rotational lines of CH + have recently been observed by 
iThi et al.l (1201 ll) in the protoplanetary disc that surrounds the 
Herbig Be star HD 100546 using the Herschel/PACS instru- 
ment. Consider ing only the excitation by nonreactive collisions, 
IThi et alJ (1201 ll) found that the strong emission lines of CH + (up 
to / = 6 - 5) arise from the rim of the outer disc and the disc sur- 
face, i.e. from media predominantly atomic with densities ~ 10 8 
cm -3 and temperatures ~ 10 2 - 10 3 K. With an effective temper - 
ature of the central star ~ 10500 K (Ivan den Ancker et a 1. 1997), 
the photodissociation is unlikely to have any effect on the dis- 
tribution of CH + among its rotational levels (see Appendix lAl). 
Conversely, the strong flux of optical photons could enhance the 
high-/ level populations at the disc surface where the CH + abun- 
dance is the highest. Moreov er, with the chemical composition 
described bv lThi et"aD d201 ll) . we expect the chemical pumping 
to dominate the nonreactive collisional processes. It would thus 
be interesting to estimate how the chemical and radiative pump- 
ing affect the understanding of the density structure of the disc 
and of the spatial emission of rotationally excited CH + . 

In preparation of future investigations, we present below an 
extension of our model to other diatomic molecules with spec- 
troscopic characteristics similar to those of CH + . 



5. Discussion: extension to other diatomic 
molecules 

5.1. Analytical considerations of the excitation timescales 

The analysis of the main excitation pathways performed in 
Appendix [B] shows that the pumping processes of CH + can be 
easily described by a set of analytical formulae. This result is in- 
herent to the X l l, + ground-state electronic configuration of CH + 
since the associated transition rules favour only one specific ex- 
citation pipeline among all others (see Fig. IB. II ). Consequently, 
the radiative pumping almost solely depends on the rotational, 
vibrational, and vibronic Einstein spontaneous emission co 
efficients, which are proportional to v 3 M; 



l 2 

'0000' 



v 3 M 2 ,„ „, and 



pumpings. In this section, we focus only on the near-infrared 
and optical/UV pumpings. 

While polyatomic molecules and other ground-state elec- 
tronic configurations will be treated in forthcoming papers (e.g. 
NH 3 , OH( 2 Il), and CO + ( 2 S + ) in priority), we thus select here 
several X'E + ground-state diatomic molecular species character- 
ized by different vibrational and vibronic Einstein coefficients: 
two light hydrides HF and HQ, and three heavy molecules 
SiO, CS, and CO. Since we aim at providing complete exci- 
tation models of these chemical compounds, we took into ac- 
count their main bound electronic configurations: the A'll states 
of CO, CS, and SiO, the state of SiO charac terized by high 
transition probabilities dElander & Smithl fl973l) . and the C'fl 
state of HC1. Conversely, we ignore the A'fl states of HF and 
HC1, which are known to correspond to repulsive configurations 
dBender & Davidson! 1968tlLedl 19991) . and the B l T. + - X'£ + sys- 
tem of HF, whose Franck-Condo n factors are found negligible 
for the v' < 10, v" = bands dAlvarino et al.lll983l) . Lastly, 
because the rotational, vibrational, and electronic structures of 
these species are clearly separated in energy, and because we 
aim to draw conclusions from their fundamental spectroscopic 
properties, we reduced their effective Hamiltonian to the three 
prime molecular constants T e (A), oj e (A), and B e (A). 

Within this approximation, the molecules considered are 
organized as follows: 



CH + SiO CS HC1 HF 

"io 6 ' " '"To 7 co To 8 



K vib = w 3 (0)M. 2 



0100 



CH + 



CS 



CO 



SiO 
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Kele = T 3 e (A)M 



A000' 



v Mi ,.. „ for the far-infrared, near-infrared, and optical/UV 



where /c v ;b and K e \ e are inversely proportional to the lifetimes of 
the excited vibrational and electronic states. According to these 
schemes and without any other considerations, CH + appears 
to be the species the least sensitive to both the near-infrared 
and optical/UV pumpings. In contrast, the vibrational transition 
dipole moment of HF and the electronic transition dipole 
moments of SiO and HC1 are so high that HF, and SiO and HC1 
are likely candidates to radiative pumpings by near-infrared 
and optical/UV photons. In the following, we compare this 
simple analytical study with the results obtained with MADEX, 
and we identify the fundamental excitation mechanisms of 
HF, HC1, SiO, CS, and CO, depending on the spectrum of 
the radiation field pervading a PDR-type medium of constant 
density «h = 10 4 cirT 3 and temperature 7k = 100 K. 

The spectroscopic properties and the nonreactive collisional 
rates of HF, HC1, SiO, CS, and CO were implemented in our 
modified version of the MADEX code, along with their respec- 
tive chemical destruction rates and photodissociation rates, both 
inferred at the border of the PDR from the Meudon PDR code. 
Except for CH + , we only considered H, H2, and He as colli- 
sion partners in the computation of the nonreactive nonelastic 
collisional rate s, thus neglecting the excitation via electron im - 
pact (see, e.g., iThummel et"allll992t lltikawa & Masonll2005l) . 
In the same way as for CH + , MADEX was then run, assum- 
ing PDR-type conditions and solar metallicities, in four differ- 
ent configurations: O considering only the excitation by non- 
reactive collisions, © considering all collisional excitation pro- 
cesses (i.e. including the chemical pumping), © considering all 
line excitation processes (i.e. including the radiative pumping 
of bound vibronic levels), and considering all radiative pro- 
cesses (i.e. including the photodissociation) with different val- 
ues of the near-infrared and optical/UV scaling factors Xw and 



10 



B. Godard et al.: A complete model of CH + rotational excitation including the radiative and chemical pumping processes 



Aopt/uv- The molecular constants, the transition dipole moments, 
the collisional rates, and the results of all models are given in 
Tables [5] [6] and [7] For comparison, we also display in Table [5] 
the results of this analysis applied to CH + . 

5.2. Species sensitive to chemical pumping 

With the chemical destruction rates considered as typical in 
PDR-type environments (see Tables [5]|7]i, we find that CH + 
is the only species, among those considered, with low-/ lev- 
els sensitive to its chemical formation. This concurs with the 
idea that chemical pumping essentially applies to short-lived 
molecules, and in particular to m olecular ions with lifetime s 
below KT 2 (10 4 /rc H ) yr (e.g. CO + . IStauber & Brudereri r2009). 
Indeed, the chemical pumping has a marginal impact on HC1, 
SiO, CS, and CO because the timescales of the nonreactive col- 
lisional excitation of the / < 5 levels are at least one order of 
magnitude lower than their respective destruction timescales. 

We note, however, that this last statement should not be gen- 
eralized to the high-./ levels. As was already stressed in Sect. 
14.51 the importance of chemical pumping relative to nonreac- 
tive collisions increases with /. For instance, Table[6]shows that 
the chemical formation takes over the excitation of the rota- 
tional levels of HF for / > 3. This suggests that this mecha- 
nism could play a dominant role for all species, especially in 
environments with a kinetic temperature close to the equiva- 
lent temperatures of the rovibrational transitions. A direct ap- 
plication could be the interpretation of the high-/ CO lines ob- 
served during the past decade in protoplanetary and planteray 
nebulae, protostellar shocks, and dense and hot PDRs with the 
KAO, SOFIA, ISO, and all the instrumen t s onboard the Herschel 
space observat ory (Justtanont et al.ll997l;lBuiarrabal et al.l2010t 
iBenedettini et alj|2012t lHabart et al.ll2010t) . If the chemical for- 
mation efficiently pumps the high-/ lines of CO, the high densi- 
ties « H > 10 6 of the CRL 618 and NGC 7027 molecular clouds 
in ferred from CO emissi on u p to the / - 22 — 21 transition 
bv lJusttanont et alJ d 19971) and Buiarrabal et all d2010h would be 
revised to lower values. The chemical pumping could similarly 
offer an alternative understanding of low-mass protostellar out- 
flows whe re the high-/ lines of CO (up to / = 38 - 37) ob- 
served by van Ke mpen et al.l d2010t) with Herschel/PACS have 
been proposed to trace the UV irradiation of the outflow cav- 
ity walls an d small-scale C-t ype shocks propagating along the 
cavity walls dVisser et alj|2012l ). 

5.3. Species sensitive to radiative pumping by near-infrared 
photons 

The comparison of models © and © shows that HF and HC1 
are the two species the most sensitive to the intensity of the 
near-infrared radiation field. In contrast, the near-infrared fluo- 
rescence is found to have almost no influence on the populations 
of the rotational levels of CH + and SiO and only a marginal ef- 
fect on those of CS and CO. When cross-compared, these be- 
haviours originate from a combination of factors: (1) the non- 
reactive collisional rates of HQ and HF (except for the 1-0 
transition) are ~ 3 to 10 times lower than those of CH + , SiO, 
CS, and CO; and (2) the lifetimes of the vibrational levels of 
HF, HC1, CO, CS, SiO, and CH + increase from HF to CH + over 
two orders of magnitude. It follows that the organization of the 
species based on their sensitivity to radiative pumping of their 
vibrational states agrees well with the classification we inferred 
from the sole analysis of the A\,ib values (see Sect. 15. lb . 



Quantitatively, the populations of the / = 1,2,3 levels of 
HF increase by a few percent to more than a factor of 10 for 
Xtm varying between 10 3 and 10 5 . A similar result is found for 
the / = 1,2,3 levels of HC1, although it requires near-inrared 
intensities ~ 10 times higher. Finally, minimum near-infrared 
radiation fields of ~ 10 6 and ~ 10 7 times that of the local ISRF 
are needed to efficiently activate the radiative pumping of CS 
and CO on one side and of SiO on the other side. 

Overall, the main tracer of near-infrared fluorescence is the 
/ = 2 - 1 transition of HF. Unfortunately, this line has, so 
far, only been observed in absorption towards Sagittarius B2, 
(lNeufeldetalJ[l997l) . In turn, the fact that the / = 1 level of 
HF is slightly enh anced from y n j T - 10 3 to 10 5 is in accordance 
with the model of I van der Tak et a il (120121) . who argued that the 
electron impact excitation is sufficient to explain the HF (1-0) 
emission line observed in the Orion Bar interclump medium with 
Herschel/HIFI. 

5.4. Species sensitive to radiative pumping of bound 
electronic states by optical and UV photons 

The predictions of MADEX concerning the excitation by opti- 
cal and UV pumping of bound electronic states are also found 
to agree well with the analysis of the K e \ e values carried out in 
Sect. 15.11 As predicted, SiO and HC1, and to a lesser extent CO, 
are sensitive to the optical and UV radiation field intensities in 
the configuration ©, while CH + and CS are not. More precisely, 
we find that the populations of the / > 3 and / > 4 levels of 
SiO and CO increase by more than a factor of 2 via the pump- 
ing of their A'n electronic states for an optical/UV radiation 
field larger than 10 3 and 10 4 (in Draine's unit) respectively. For 
Xopi/uv - 10 4 , we find that this process even induces a popula- 
tion inversion of the / = and / = 1 rotational levels of SiO. 
This remarkable behaviour contrasts with that of CH + and CS, 
whose rotational level populations are influenced by optical and 
UV pumping only for a radiation field 10 s times larger than the 
local ISRF. 

The applications of the radiative pumping of bound elec- 
tronic states of SiO extend to a variety of astrono mical sources- 
Obser ved in the Orion Bar and S 140 PDRs by ISchilke et all 
(120011) . the high-/ (up to / = 5 - 4) and low-/ (/ = 2 - 1 and 
/ = 3 — 2) rotational lines of SiO have been thought to trace 
clumps with densities of 2 x 10 6 cirT 3 and molecular clouds at 
lower densities (down to a few 10 4 cm -3 ), respectively. With 
the predictions of Table [5] we estimate, on the contrary, that the 
entire SiO emission originates from the low-density component 
and that the high-/ lines are confined to the border of the PDR 
where the radiative pumping is maximal. In addition to repro- 
ducing the SiO rotational diagram, this last result is in concor- 
dance with the extended nature of the SiO (2-1) emission - 
compared to that of SiO (5 - 4) - observed in the Orion Bar 
PDR. 

Similar conclusions may be drawn for the high-velocity 
molecular outflows of class protostars, where the high-/ (up 
to 11 - 10) rotational lines of SiO associated to the molecular 
bullets distant from the central star are detected with intensi- 
ties considerably lowe r than those measu r ed in the bullets close 
to the driving source dNisini et al.l l2002[. [2007). Originally at- 
tributed to density gradients (Nisi ni et al.ll2007l) . this behaviour 
may be explained with the dilution of the strong FUV radiation 
field emitted by the young stellar object along the outflow cavity 
(e.g. lBruderer et al.ll2009h . 

Finally, the radiative pumping by optical and UV photons 
could play a major role on the production of the SiO masers 
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observed in the region between the photosph ere and the dust- 
growth zone of evolve d stars envelopes (e.g. iGrav et al.llT995t 
ICernicharo et aljfl993l and references therein). With the transi- 
tion dipole moments of SiO given in Table [5] this process is ex- 
pe cted to drive the e xcitation of the v > 2 levels, in contradiction 
to Buiarrabal ( 1994), who proposed direct infrared pumping of 
the Av = 1 transitions as a global inversion process. In a sub- 
sequent paper, we will confront the predictions of our complete 
excitation model to the forthcoming survey of SiO maser emis- 
sion (up to v = 6) performed on a sample of 50 O-rich evolved 
stars (in preparation). 

5.5. Species sensitive to their photodissociation 

Finally, Tables [5] |6] and [7] show that the photodissociation pro- 
cesses have a substantial effect on the rotational diagrams of all 
species considered for^ opt /uv varying between 10 2 and 10 4 . In 
particular, the photodissociation, whose rates are derived for an 
optical/UV flux scaling with the standard ISRF, is found to com- 
pete with the bound-bound absorption of optical/UV photons by 
SiO, CO, and CH + , and even to dominate the bound-bound ab- 
sorption of optical/UV photons by HC1 and CS. 

The species the most sensitive to this mechanism are HF, 
HC1, and to a lesser extent CS. For instance, we find that the 
photodissociation drives the populations of the J > 4 levels of 
HF and HC1 and of the J > 3 levels of CS, providing that the gas 
is illuminated by a UV radiation field 10 2 and 10 3 times that of 
the local ISRF. Moreover, since the photodissociation acts as an 
additional chemical destruction rate, its impact on the level pop- 
ulations is similar to the chemical pumping and thus increases 
with J. 

These results suggest that the photodissociation may influ- 
ence the rotational diagrams of all the species singled out in 
Sect. 15.21 15.31 and 15.41 We note, however, that the importance 
of the process relative to the chemical pumping and the radiative 
pumping of bound states strongly depends on the chemical com- 
position of the gas and the spectrum of the optical/UV field. We 
therefore recommend to use models that compute the coupling 
between the radiative transfer and the chemical evolution of the 
gas to accurately describe the excitation of molecules detected 
in bright PDRs. 



6. Conclusions and perspectives 

We have performed a theoritical investigation of the main exci- 
tation and de-excitation mechanisms of the methylidyne cation 
in several astrophysical environments. To do so, we built a new 
version of the MADEX radiative transfer model that not only 
includes all rotational, rovibrational, and rovibronic transitions 
of CH + , but also takes into account the excitation induced by 
reactive and nonreactive collisions. The model was used to es- 
timate the optical and near-infrared radiation fields required to 
perturb the steady-state rotational level populations of CH + , and 
also to provide valuable predictions of the intensities of its far- 
infrared emission lines. The results of MADEX were systemati- 
cally compared with an extensive analytical study of the pump- 
ing processes of CH + , which we finally used to extend our pre- 
dictions to five other ground-state diatomic molecules, HF, 
HC1, SiO, CS, and CO. 

A remarkable agreement was found, throughout this paper, 
between the results of MADEX and the analytical description, 
thus giving access to the dependence of the excitation processes 
on the main physical conditions of the gas («h, Tk, and I v ) and on 



the fundamental spectroscopic properties of the molecule con- 
sidered. The key parameters that govern the strength of the ra- 
diative pumping of the vibrational and bound 1 n electronic states 
of a given species are /c v ib = w 3 (0)M„ 1()0 and K de = r 3 (l)M 2 0(K) . 
It follows that the far-infrared and submillimetric lines of HF 
and HC1 (SiO, HC1 and CO) strongly depend on the intensity 
of the near-infrared (optical/UV) radiation field in astrophysical 
environments submitted to a photon flux at least 10 3 x («h/10 4 ) 
and 10 4 x («h/10 4 ) times the value observed in the solar neigh- 
bourhood. Such environments include bright PDRs where the 
photodissociation is expected to compete with the line pro- 
cesses in the rotational excitation of the molecules. The species 
the most sensitive to the photodissociation are HF and HC1, 
whose rotational diagrams are affected by a UV flux as low as 
10 2 X (« H /10 4 ) times that of the local ISRF. In constrast, the dis- 
tribution of populations among the rotational levels of CH + is 
weakly sensitive to both the radiative pumping of bound states 
and the photodissociation, and is found to strongly depend on 
the CH + chemical formation. 

As a natural consequence, the intensities of the pure rota- 
tional lines of CH + observed towards the Orion Bar and the NGC 
7027 PDRs are reproduced with MADEX assuming molecular 
clouds with proton densities of 5 x 10 4 and 2 x 10 5 cirT 3 , respec- 
tively. Our estimate of the density of the NGC 7027 circumstellar 
PDR is one to two orders of magnitude below the value inferred 
from traditional excitation models. In addition, we found that 
the CH + emission originates from the interclump medium of the 
Orion Bar PDR, which contradicts the spatial correlation found 
between the CH + (3 - 2) emission line and the high-density trac- 
ers such as the rotational A-doublets of OH and high-/ rotational 
transitions of CO. Consequently, all our results favour molecular 
clouds with densities significantly lower than previously estab- 
lished. 

It is foreseeable that similar investigations will have to be 
carried out for all molecules detected in interstellar and circum- 
stellar molecular clouds. In forthcoming papers we will address 
in priority the following problems. 

1. The radiative pumping of the A l Tl and electronic lev- 
els of SiO may have a dominant role in the production of 
SiO masers in the envelopes of evolved stars. We intend to 
perform a statistical analysis of the inversion process based 
on the comparison of the predictions of MADEX with the 
forthcoming observations of the SiO masers in the high vi- 
brational sates (up to u = 6) in a sample of 50 O-rich evolved 
stars covering a wide range of mass loss rates. 

2. Because the efficiency of the radiative pumping highly de- 
pends on the electronic configuration of the considered 
species, we plan to extend the predictions of MADEX to 3 E + 
and 2 n ground-state diatomic molecules. Our goal is to study 
two molecules that are believed to trace media with high den- 
sity and temperature, CO + and OH. 
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Appendix A: Analysis of the dominant de-excitation 
processes of CH + 

The concept of critical density - density at which the nonre- 
active collisional de-excitation rate equals the radiative decay 
rate - provides valuable insight into the mechanism that con- 
trols the population of each level of a molecule depending on 
the physical conditions of the gas. However, this concept is in- 
sufficient if state-to-state chemistry is taken into account, or if 
the wavelengths of the radiative de-excitations of the level con- 
sidered span a wide range of the electromagnetic spectrum, as is 
the case if rovibrational or rovibronic transtitions are considered. 
We therefore address the following question: in this general case, 
can we derive a set of critical parameters that will give us insight 
into what the steady-sate level populations should be depending 
on the conditions of the gas ? 

We consider a species S whose levels are characterized by 
three quantum numbers only, the total orbital angular momen- 
tum A, the vibrational momentum v, and the total angular mo- 
mentum /. Taking into account all radiative and collisional pro- 
cesses, the density of atoms «, in the level i (with a degeneracy 




le+05 



Figure A.l. Dominant de-excitation process of CH + J — 1 as a 
function of the density «h, and the equivalent black body tem- 
peratures r v jb and r e i e , for T Iot = 2.7 K. In the blue, purple, red, 
and brown regions, the de-excitation of CH + J — 1 is dominated 
by pure rotational, rovibrational, and rovibronic radiative transi- 
tions and by the photodissociation. In the empty space above, the 
de-excitation occurs through reactive and nonreactive collisions 
depending on and the chemical composition of the gas (see 
main text). 



gi) verifies, at steady-state, 

F c I j «0 

n, = ; jr-^ 

D C j j S' J 

(A.l) 

and 

«(S) = ^ "i (A- 2 ) 

In these formulae, nr>, rif a , and rif h are the densities of reagents 
that chemically destroy and form the level i at rates kP and k? , 
nc and kfj are the abundances of the nonreactive collision part- 
ners and the corresponding collisional rates, Ay and (n y . v ) are 
the Einstein spontaneous de-excitation coefficients and the di- 
rectionally averaged photon occupation numbei0of the radiation 
field at the frequency v of the i - j transition, and cr 1 ! are the pho- 
todissociation cross-sections of the level i. 

To estimate the importance of the different de-excitation pro- 
cesses, we deliberately decomposed the radiative term of Eq. lA.ll 
into three parts that include all transitions: (1) between two lev- 
els in the same electronic and vibrational state (A, = Aj and 
Vj = vj), (2) between two levels in the same electronic state 
and different vibrational states (A,- = Aj and u, + vj), and (3) 
between two levels in different electronic states (A,- + Aj). We 
then regrouped all processes involving high-energy photons, i.e. 
the bound-bound rovibronic transitions and the photodissocia- 
tion. In addition, we assumed that the radiation field intensities 
at the corresponding frequencies are described by three black 
bodies with equivalent temperatures T rot , T v n,, and T e \ e respec- 



4 ( n y,v) is related to the radiative energy density u v (in erg cirT 3 Hz~') 
and the radiative specific intensity I v (in erg cur 2 Hz 1 sterad -1 s _1 ) 
through : (n yy ) = /„c 2 /2/iv 3 and u v = 4nl v /c 
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tively. Re-arranging Eq. IA.ll we obtain 



(A.7) 



Tj"F a n F b *f 




1 + • 



"H crij 



Z =u >>j(Aji + (ny,v)Aji + f.<Py, v )A is ) 



2^. (Aij + iny^Atj+^iny^Aj,) 



(A.3) 



X n } (A ji + {n„)Aji + j- (n y , v )Aij) 



1 + 



Zerij 



1 + • 



1 + • 



E ji + <« r ,v>Aji + f. (n„)Aij) 

Z (Ay + <n r ,y>Ay + Si<n r ,v>A^) ' 

where each term tends to drag the level populations towards dif- 
ferent Boltzmann distributions at temperatures TV, 7k, T I0t , T v n,, 
and r e i e . In the previous equation we define five independant 
parameters, the critical total hydrogen density «h cri,/ (in cm" 3 ), 
and four dimensionless critical ratios, X c „ j which depend on the 
abundances relative to H of the reactive and nonreactive colli- 
sional partners [D] and [C], and Fori,!, Z cr y, and S C ry, which de- 
pend on the entire spectrum of the radiation field: 



«Hc 



£(Ay + <« r ,v)Ay + g Mn %v )Afi) + J ^[WWv 
j 

ZlD]k? + ZZlC]kf 

D C j J 



c j 



* crU = 2[D]*f 



(A.4) 



(A.5) 



y ■ = 

■*- cx\ .i — 



2 (Ay + (n 7iV >Ay + ff <n r , v >A ji ) + 

2 (Ay + ( % v)A,i + |<» r ,v)A i7 ) + / %crj(v)l v dv 



2 (Ay + <n r>v >Ay + ^<n y , v )A J i) 

■// A j = A,- 



(A.6) 



7 - 



2 (Ay + (« r , v )Ay + |<n y ,v)A i7 ) + / f y crj(v)l v dv 
2 (Ay + (n r , v )A, 7 + ff<n T , v >A^) 

Jf A,=A, 



5" cri.i — 



2 (Ay + (« r ,y>Ay + f-(n %v )Aji) 



(A.8) 



These quantities prove useful to explore the parameter domain 
since they indicate the main de-excitation mechanism of each 
level depending on the condition of the gas: through reactive col- 
lisions for «h > «h cri,; and X or y < 1, nonreactive collisions for 
«h > «h en,; and X a -ij > 1, rotational radiative transitions for 
«h < «Hcri,! and Ycnj < 1, vibrational radiative transitions for 
«h < «Hcri,(, lay > 1 and Z cr j,, < 1, electronic radiative tran- 
sitions for «h < «Hcri,i, Y cr ij > 1, Z cl y > 1, and S cr ij > 1, 
and photodissociation for «h < «Hcii,i, T cr y > l,Z CI y > 1, and 

S cri,;' < 1- 

When applied with the spectroscopic, collisional and chem- 
ical properties of CH + , the critical parameter analysis gives the 
following results: 

• given their short lifetimes, the de-excitation of the vibra- 
tional and electronic levels of CH + always occurs through 
spontaneous radiative emission of near-/mid-infrared and op- 
tical photons. 

• The critical ratios X a [j < 1 for all rotational levels J > 1 
as long as T K < 710 x (2.5 - 3.0[H 2 ])" L9 K and T K < 
1 1 00 x ( [e~] / 1 (T 4 )~ 1 K. The other critical parameters mainly 
depend on the density of the gas and the intensity of the ra- 
diation field. 



As an example, we therefore display in Fig. lA.ll the limits set 
by the critical parameters analysis for CH + in the level A,v,J = 
0,0,1, assuming [H 2 ] = 0.5, [He] = 0.1, [e~] = 10~ 4 , T K = 
100 K, 7V = To = Tk, and T mt = 2.7 K. The edges at «h ~ 
10 7 cm" 3 , r vib ~ 630 K, T ele ~ 1740 K, and r e i e ~ 14400 K 
point to the physical conditions for which several de-excitation 
mechanisms have to be taken into account to accurately compute 
the level populations. 

According to Fig. IA.ll a large UV and optical radiation field 
seems to be required for the photodissociation to compete with 
the bound-bound rovibronic absorption. However, we emphasize 
that this result strongly depends on the functional form used to 
described the radiation field, i.e. a black body with an effective 
temperature T e \ e in the present case. In the standard unshielded 
interstellar radiation field, the photodissociation and the bound- 
bound rovibronic absorption are found to have comparable rates 
and are therefore expected to play similar roles in distributing 
CH + among its rotational levels. 

Appendix B: Analysis of the dominant excitation 
processes of CH + 

Fig. IB. II shows the excitation and de-excitation pathways that 
drive the populations of the first five rotational levels of CH + 
computed with MADEX and neglecting the photodissociation. 
Five extreme models were selected: (a) 7k = 100 K, x&r = 1, 
Xmt = 1, and^ p t = 1, (b) T K = 1000 K,^ fc = l,Xmc = 1, and 
Xov x = 1, (c) T K = 100 K, X nr = 10 6 ,*nir = 1, and^ opt = 1, (d) 
T K = 100 K,Xfk = hXm = 10 10 , and^ opt = 1, (e) T K = 100 
K,*fir = hXmr = 1, and xopt = 10 6 , assuming 7V(CH+) = 10 13 
cm -2 , v exp = 1 km s _1 , and «h = 10 4 ctrT 3 . This figure is simpli- 
fied: for each level, only the processes that together contribute to 
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(a) 




(b) 



reagents 
C + + H 2 
CHj + y 



products 
CHj + H 



(C) 



(d) 




v=0 v=l 




10 13 cm- 2 , 



^exp 



1 km s 



Figure B.l. Main excitation and de-excitation routes computed with MADEX, assuming N(CH + ) 

« H = 10 4 cm- 3 , and (a) T K = 100 K, Xsi = = h and^ opt = 1, (b) T K = 1000 K,^ fir = l, X w = h and^ opt = 1, (c) T K = 100 

K, Xfir = 10 6 ,x™ = 1, and^ opt = 1, (d) T K = 100 K,^ fir = \, X w = 10 10 , and* opt = 1, (e) T K = 100 K,* fir = l,* nir = 1, and 
Xopt = 10 6 . The photodissociation process was intentionally neglected to focus on the (de)excitation of the bound-bound rovibronic 
transitions. The pure rotational, rovibrational, and rovibronic transitions are displayed in blue, purple, and red. Nonreactive and 
reactive collisional processes are displayed in green and yellow. In each case, the levels that do not contribute to the excitation 
of the pure rotational levels are not displayed (e.g. A,v,J = 0, 1, in the case (d)). For more clarity the rotational emission lines 
from 7 = 10 to 7 = 5 are merged into a single dotted arrow in case (b), and all vibronic A, 7 = 1,1 <— 0, absorption lines and 
A, 7 = 1,1 — * 0, 2 emission lines are merged into two single arrows in case (e). 



more than 70 percent of the total destruction and formation rates 
are displayed. 

In all cases the de-excitation of the rotational levels occurs 
via pure rotational transitions, as expected from the analytical 
treatment presented in appendix [A] Conversely, four different 
excitation regimes are revealed: (a) by reactive collisions for all 
7 > 1 levels, (b) by chemical pumping of a high-energy level 
(e.g. 7=10) followed by the cascade through the intermediate 7 
levels, (c) by absorption of far-infrared photons by the 7-1 level, 
(d) by near-infrared pumping of the A, v = 0, 1 levels, and (e) 
by optical pumping of the A, 7 = 1,1 levels followed by the cas- 
cade through the vibrational levels of the X l I. + state. Using Eq. 
lA.ll and taking into account all routes displayed in Fig. IB. II we 
derive at the steady-state equilibrium 



»007 

n(CH+) 



= «H" 



([H]^ + [H 2 ]^ 2 + [ e -]^) 



A)oyoo(y-i) 
(27 + 1) exp(-£ooj/r K ) 
2(2/ + 1) exp(-£,/r K ) 



for 7 > 1 in case (a), 

([H]£g + [H 2 ]£g 2 + [ e -]^) 



(B.l) 



"007 

n(CH+) 



= «H" 

^00700(7-1) 

max [(2/ + l)exp( -.E, / TV)] 

!>7 

2(2/+l)exp(-£,/r F ) 



(B.2) 



for 7 > 3 in case (b), 



«007 =«000 (27+ 1) Xfo f|(«r,v a -i> 



1=1 



for 7 > 2 in case (c), 



«007 =«ooo X, 
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\~[ ^00(2; 



;-l)00(2/-2) 



V (=1 



(V-D/2 

n 



4i+l («y,v)Aoi(2i)00(2i-l) 



U 4M 1 + 



>ioi(2000(2(+]) / 



, («y,vM011000 

1 _^ ^01 1000 

A011002 / 



for odd 7-values in case (d), 
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for even /-values in case (d), and 

Xopt 

«oo./ =»ooo -; x 

Z„ \ A ^li/10(7-2)2 1 
s{n 7 , v )A lvim) —- x 
v v 2jv> AiuWi/O ) K ub ) 

I 1 . \ 



111, ^0(J-i+l)(i-l)0(J-.)(.-2) 

;_-> 1 -T — T 

V ' J «0(/-j'+l)(i-l)0(/-i)i / 

for J > 2 in case (e). In all these formulae, (« 7>v > corresponds 
to the directionally averaged photon occupation number of the 
standard interstellar radiation field at the frequency v of the ra- 
diative transitions. 



B. Godard et al.: A complete model of CH + rotational excitation including the radiative and chemical pumping processes 
Table 5. Predictions of MADEX applied to CH + and SiO. 



"CH + 



"SiO 



T e (cm" 1 ) 
2.412 (+04) 



u> e (cm ) 

2.858 (+03) 
1.864 (+03) 



B e (cm" 1 ) 

1.418 (+01) 
1.189 (+01) 



T e (cm" 1 ) 



4.284 (+04) 



u) e (cm [ ) 

1.242 (+03) 
8.528 (+02) 



B e (cm" 1 ) 

7.242 (-01) 
6.307 (-01) 

















M A > 


j'A" 


v n (eag) 












„" (eao) 








v"\u' 







1 




2 




3 


4 







1 


2 


3 


4 









2.2 (- 


-1) 


1.6 (- 


-1) 


1.0 (- 


-i) 


6.5 (-2) 


4.1 (- 


-2) 


2.0 (+0) 


2.7 (+0) 


2.6 (+0) 


2.0 (+0) 


1.3 (+0) 


< 


c 


1 


1.1 (- 


-1) 


9.1 (- 


-2) 


1.5 (- 


-i) 


1.3 (-1) 


1.1 (- 


-1) 


2.7 (+0) 


1.9 (+0) 


1.6 (+0) 


1.5 (+0) 


2.1 (+0) 


T 


T 


2 


4.6 (- 


-2) 


1.1 (- 


-1) 


1.6 (- 


-2) 


7.0 (-2) 


1.0 (- 


-1) 


2.8 (+0) 


1.4 (-1) 


1.8 (+0) 


1.5 (+0) 


6.6 (-2) 


< 




3 


1.6 (- 


-2) 


7.2 (- 


-2) 


6.5 (- 


-2) 


6.7 (-2) 


1.5 (- 


-2) 


2.4 (+0) 


1.3 (+0) 


1.7 (+0) 


3.4 (-1) 


1.6 (+0) 




4 


5.1 (- 


-3) 


3.5 (- 


-2) 


7.6 (- 


-2) 


2.1 (-3) 


5.2 (- 


-2) 


1.8 (+0) 


2.1 (+0) 


4.8 (-1) 


1.6 (+0) 


9.8 (-1) 
















M x , 


j'A" 


v" (ea ) 










MA'v'A' 


v" (ea Q ) 








u"\u' 







1 




2 




3 


4 







1 


2 


3 


4 









6.6 (- 


-1) 


6.3 (- 


-3) 


6.6 (- 


-5) 


6.5 (-7) 


6.5 (- 


-9) 


2.0 (+0) 


4.1 (-2) 


2.3 (-3) 


1.6 (-4) 


1.9 (-5) 


< 


c 


1 






6.6 (- 


-1) 


9.6 (- 


-3) 


7.7 (-4) 


2.0 (- 


-5) 




6.1 (+0) 


5.7 (-2) 


4.4 (-3) 


3.2 (-4) 


T 


T 


2 










6.7 (- 


-1) 


1.3 (-2) 


2.9 (- 


-3) 






1.2 (+1) 


6.9 (-2) 


6.2 (-3) 


< 


c 


3 














6.7 (-1) 


1.8 (- 


-2) 








1.8 (+1) 


7.8 (-2) 




4 
















6.8 (- 


-1) 










2.0 (+1) 



j"\r 
o 
i 

2 
3 
4 
5 





1.6 (-10) 
8.0 (-11) 
7.0 (-12) 

1.4 (-12) 

2.5 (-14) 



1 



£[C]^„ (cm 3 s- 
2 3 



8.1 (-11) 5.3 (-11) 1.1 (-11) 8.4 (-12) 9.1 (-13) 



6.6 (-11) 

2.4 (-11) 

9.5 (-13) 
2.5 (-13) 



8.8 (-11) 

2.5 (-11) 
9.2 (-12) 
9.8 (-14) 



7.5 (-11) 
5.9 (-11) 

9.6 (-12) 
3.1 (-12) 



1.2 (-11) 
8.4 (-11) 
3.7 (-11) 

4.0 (-12) 



1.9 (-11) 
5.4 (-12) 
7.1 (-11) 
2.4 (-11) 





8.4 (-11) 
5.0 (-11) 

4.5 (-11) 
8.5 (-12) 
1.7 (-11) 



1 

2.8 (-11) 

7.8 (-11) 
3.4 (-11) 

3.4 (-11) 

5.5 (-12) 



£[C]fc£ 7 „ (em's- 1 ) 
2 3 



1.1 (-11) 
4.9 (-11) 

7.1 (-11) 
2.9 (-11) 
2.8 (-11) 



7.3 (-12) 
1.6 (-11) 

5.4 (-11) 

6.6 (-11) 
2.6 (-11) 



1.2 (-12) 

1.4 (-11) 
1.9 (-11) 

5.5 (-11) 

6.1 (-11) 



2.1 (-12) 
2.0 (-12) 
1.7 (-11) 
2.0 (-11) 
5.6 (-11) 



2 [D]k D (cm 3 s- 1 ) 
6.24 (-10) 



/ %<rv(v)l v dv (s- 1 ) 
3.30 (-10) 



lZ[D]k D (cm 3 s- 1 ) 
9.35 (-14) 



/ %o-?(v)I v dv (s- 1 ) 
1.60 (-09) 



./ 
O 
© 







1 



2 3 







1 



«j/n(SiO) 
2 3 



10 3 

£ io 5 

* io 7 

I io; 

a 10 3 

% 10 4 

i 10 

X 10 4 









1.0 (+00) 1.4 (-03) 1.5 (-05) 3.8 (-07) 2.6 (-08) 3.2 (-10 

1.0 (+00) 4.3 (-03) 5.9 (-05) 4.8 (-06) 4.5 (-07) 3.4 (-08 

1.0 (+00) 4.3 (-03) 5.9 (-05) 4.8 (-06) 4.5 (-07) 3.4 (-08 

1.0 (+00) 4.3 (-03) 5.9 (-05) 4.8 (-06) 4.5 (-07) 3.4 (-08 

1.0 (+00) 4.6 (-03) 7.0 (-05) 4.8 (-06) 4.5 (-07) 3.4 (-08 

1.0 (+00) 4.3 (-03) 5.9 (-05) 4.8 (-06) 4.5 (-07) 3.4 (-08 

1.0 (+00) 4.3 (-03) 6.1 (-05) 4.9 (-06) 4.5 (-07) 3.4 (-08 

9.9 (-01) 5.0 (-03) 7.2 (-05) 5.4 (-06) 4.8 (-07) 3.6 (-08 

1.0 (+00) 4.3 (-03) 6.0 (-05) 4.8 (-06) 4.5 (-07) 3.4 (-08 

1.0 (+00) 4.5 (-03) 6.3 (-05) 5.1 (-06) 4.7 (-07) 3.6 (-08 

9.9 (-01) 6.4 (-03) 9.9 (-05) 8.0 (-06) 7.2 (-07) 5.5 (-08 



2.6 (-01) 4.5 (-01) 2.3 (-01) 4.2 (-02) 3.5 (-03) 3.1 (-04) 
2.6 (-01) 4.5 (-01) 2.3 (-01) 4.2 (-02) 3.5 (-03) 3.1 (-04) 



2.6 (-01) 4.5 (-01) 2.3 (-01) 4.2 (-02) 

2.6 (-01) 4.5 (-01) 2.3 (-01) 4.2 (-02) 

2.1 (-01) 4.2 (-01) 2.8 (-01) 7.1 (-02) 

2.5 (-01) 4.5 (-01) 2.3 (-01) 4.5 (-02) 

2.1 (-01) 4.3 (-01) 2.6 (-01) 7.1 (-02) 

9.3 (-02) 3.0 (-01) 3.1 (-01) 1.9 (-01) 

2.5 (-01) 4.5 (-01) 2.3 (-01) 4.8 (-02) 
1.9 (-01) 4.1 (-01) 2.8 (-01) 9.0 (-02) 

6.6 (-02) 2.5 (-01) 2.9 (-01) 2.2 (-01) 



3.5 (-03) 3.1 (-04) 

3.6 (-03) 3.1 (-04) 

7.2 (-03) 5.4 (-04) 

4.0 (-03) 3.7 (-04) 

9.3 (-03) 1.1 (-03) 

7.1 (-02) 1.8 (-02) 

4.5 (-03) 4.8 (-04) 

1.5 (-02) 2.4 (-03) 

1.1 (-01) 3.8 (-02) 



From top to bottom: molecular constants of the X l Y, + and A'll electronic states, transition dipole moments of the A l Yl — X l 2. + and X l ~L + -X [ l. + 
band systems, nonreactive collisional rates, chemical destruction rates and photodissociation rates, and populations of the rotational levels 
predicted with MADEX for J < 5. The collisional rates and the predictions of the level populations are computed for T D = 7p = T K = 100 K, 
«h = 10 4 crrr 3 , [Ho] = 0.5, and [e~] = 10~ 4 . For both molecules, MADEX was run in the four configurations O, ©, ©, and (see main text). 
In cases © and 0, the values of the radiation field scaling factors Xmi an d^ op t/uv are 1 unless specified otherwise. Numbers in parenthesis are 

power of ten. 

data from lHakalla et alj <2006l); ICheng etall ( 120071): iLim et alj dl999 T> ; iHammami et al.l d2008h : iTurpin et al.| d2010h: fc da t a from 
Hedelund & Lambed dl972h : iBeer et alJ dl974h: Barrow & Stond il975h: lElander & SmifhMl973h: iNaidul d!98lh : iDrira et~aH d!997l) : 
;[Lii 



Tipping & Chackeriarl i fl981h : lLiszt & SmitH dl972h ; lDrira et alj dl998l) : lPark et al.Ul993l) ; lDavou & Balancal d2006) 
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B. Godard et al.: A complete model of CH + rotational excitation including the radiative and chemical pumping processes 
Table 6. Same as Table|5]for HF and CS. 



"HF 



"CS 



A 


1 



T e (crrT 1 ) 



oj e (cm ') 
4.138 (+03) 



B e (cm 4 ) 
2.056 (+01) 



T e (cm" 1 ) 
3.890 (+04) 



u) e (cm ) 

1.285 (+03) 
1.073 (+03) 



B e (cm" 1 ) 

8.200 (-01) 
7.800 (-01) 



v"\v' 


1 

2 
3 
4 



v"\u' 


1 

2 
3 
4 



2 3 



A/ A ^' A «„« (ea ) 








1 




2 




3 




4 




2.0 (- 


-1) 


9.2 (- 


-2) 


2.9 (- 


-2) 


6.6 (- 


-3) 






6.4 (- 


-2) 


1.6 (- 


-1) 


1.2 (- 


-1) 


5.2 (- 


-2) 


3.6 (- 


-2) 


1.9 (- 


-2) 


8.3 (- 


-2) 


1.2 (- 


-1) 


1.3 (- 


-1) 


2.2 (- 


-1) 


4.8 (- 


-3) 


3.3 (- 


-2) 


9.1 (- 


-2) 


8.3 (- 


-2) 


4.7 (- 


-1) 


7.8 (- 


-3) 


6.2 (- 


-2) 


2.4 (- 


-1) 


4.2 (- 


-1) 


1.3 (- 


-1) 





7.2 (- 



M A ' u 'A"v" (ea ) 
2 3 



1) 3.6 (-2) 
7.3 (-1) 



4.5 (-3) 
4.7 (-2) 
7.5 (-1) 



5.6 (-4) 
7.5 (-3) 
5.3 (-2) 

7.7 (-1) 



8.2 (- 
1.1 (- 
1.0 (- 
5.6 (- 
7.8 (- 





7.7 (- 



2 3 



1) 6.2 (-2) 
7.6 (-1) 



3.6 (-3) 
8.8 (-2) 
7.4 (-1) 



2.4 (-4) 
6.3 (-2) 
3.9 (-2) 
7.3 (-1) 



2.0 (-5) 

2.0 (-3) 
1.2 (-2) 
3.9 (-2) 

7.1 (-1) 



J"\J' 


1 

2 
3 
4 
5 



1.2 (-10) 

2.7 (-12) 
1.5 (-13) 

3.8 (-15) 

1.3 (-16) 



1 

7.0 (-11) 

2.5 (-11) 

1.6 (-13) 

3.7 (-14) 
7.6 (-17) 



Z[C]^„ 
2 

3.3 (-12) 
4.9 (-11) 

5.4 (-12) 
2.3 (-14) 
4.6 (-15) 



(cm 3 s ') 
3 

7.9 (-13) 
1.4 (-12) 

2.4 (-11) 

1.5 (-12) 
1.4 (-13) 



4 

1.8 (-13) 
2.8 (-12) 
8.7 (-13) 
1.3 (-11) 

5.1 (-13) 



9.7 (-14) 
9.6 (-14) 
2.9 (-12) 
2.1 (-11) 
8.6 (-12) 



4.9 (-11) 
8.6 (-11) 
1.9 (-11) 
2.1 (-11) 
5.8 (-12) 



1 

1.7 (-11) 

4.0 (-11) 
5.9 (-11) 

1.3 (-11) 

1.4 (-11) 



Y,[C]k c JtJ „ (cm 3 s- 1 ) 
2 3 



1.9 (-11) 
2.5 (-11) 

3.4 (-11) 
4.8 (-11) 
1.1 (-11) 



3.1 (-12) 
2.9 (-11) 
2.6 (-11) 

3.3 (-11) 
4.0 (-11) 



4 

2.9 (-12) 
5.3 (-12) 
3.2 (-11) 
2.8 (-11) 

2.8 (-11) 



7.5 (-13) 

5.4 (-12) 

6.5 (-12) 
3.2 (-11) 

2.6 (-11) 



2 [D]k D (cm 3 s" 1 ) 
8.42 (-13) 



/ %o-v(v)I v dv (s- 1 ) 
1.17 (-10) 



Z[D]k D (cm 3 s" 1 ) 
1.16 (-14) 



/ fy(V)l v dv (s- 1 ) 
9.80 (-10) 



»y/»(HF) 



./ 
O 
© 





1.0 (+00) 

1.0 (+00) 

1.0 (+00) 
1.0 (+00) 
9.9 (-01) 



1 

2.0 (-04) 

2.0 (-04) 

2.0 (-04) 
2.8 (-04) 
8.2 (-03) 



10 2 

10 3 

10 4 



1.0 (+00) 2.0 (-04) 
1.0 (+00) 2.1 (-04) 
1.0 (+00) 3.4 (-04) 



2 

1.2 (-07) 

1.3 (-07) 

1.5 (-07) 
2.1 (-06) 
2.0 (-04) 



1.4 (-07) 
2.6 (-07) 

1.5 (-06) 



3 

1.7 (-09) 

2.3 (-09) 

2.4 (-09) 
1.1 (-08) 
1.1 (-06) 



3.0 (-09) 
9.9 (-09) 
7.8 (-08) 



4 

1.8 (-11) 
4.3 (-11) 

4.3 (-11) 

5.9 (-11) 
5.3 (-09) 







1 



nj/n(CS) 
2 3 



2.9 (-13) 

1.0 (-12) 

1.0 (-12) 
1.0 (-12) 
2.5 (-11) 



2.5 (-01) 4.8 (-01) 2.2 (-01) 3.9 (-02) 5.3 (-03) 1.2 (-03) 
2.5 (-01) 4.8 (-01) 2.2 (-01) 3.9 (-02) 5.3 (-03) 1.2 (-03) 



7.6 (-11) 2.0 (-12) 
3.8 (-10) 1.1 (-11) 
3.4 (-09) 9.9 (-11) 



2.5 (-01) 4.8 (-01) 

2.5 (-01) 4.8 (-01) 

1.3 (-01) 3.9 (-01) 

2.5 (-01) 4.8 (-01) 

2.5 (-01) 4.8 (-01) 

2.4 (-01) 4.7 (-01) 

2.5 (-01) 4.7 (-01) 
2.1 (-01) 4.5 (-01) 
6.8 (-02) 2.9 (-01) 



2.2 (-01) 
2.2 (-01) 
3.2 (-01) 



3.9 (-02) 

4.1 (-02) 

1.2 (-01) 



5.3 (-03) 1.2 (-03) 

5.4 (-03) 1.2 (-03) 

2.5 (-02) 4.0 (-03) 



2.2 (-01) 3.9 (-02) 

2.2 (-01) 4.0 (-02) 

2.3 (-01) 4.3 (-02) 

2.2 (-01) 4.2 (-02) 

2.5 (-01) 6.2 (-02) 

3.2 (-01) 1.8 (-01) 



5.3 (-03) 
5.3 (-03) 
5.8 (-03) 

6.0 (-03) 
1.3 (-02) 
7.2 (-02) 



1.2 (-03) 
1.2 (-03) 
1.2 (-03) 

1.5 (-03) 
4.2 (-03) 
3.0 (-02) 



' dat a fromlWebb & Narahari Rad dl968l): Uohns & Barrow! _ . 

from lLovas & Krupenid dl974l) : IWinnewisser & Cookl d 19681) ; IB otschwina & SebalJ jl985h : lTurner et alj 



1959h:lMeredithl dl973l)^lReese et al] d2005h : 



data 



Guillon & Stoecklinl d2012h 
1 1992h ; lLique & S pielfiedel d2007h 
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B. Godard et al.: A complete model of CH + rotational excitation including the radiative and chemical pumping processes 

Table 7. Same as Table[5]for HC1 and CO. In contrast to other molecules, the A = 1 state of HQ corresponds to its C'fl (and not 
A 1 !!) electronic configuration. 



"HC1 



"CO 



A 


1 



T e (cm" 1 ) 
7.758 (+04) 



u) e (cm ') 

2.990 (+03) 
2.684 (+03) 



B e (cm" 1 ) 

1.059 (+01) 
9.330 (+00) 



T e (cm- 1 ) 
6.508 (+04) 



to e (cm l ) 

2.170 (+03) 
1.518 (+03) 



B e (cm" 1 ) 

1.931 (+00) 
1.612 (+00) 



v"\v' 


1 

2 
3 
4 



v"\u' 


1 

2 
3 
4 





8.6 (- 

5.8 (- 

4.2 (- 

1.9 (- 

1.3 (- 





4.3 (- 



Mavav (ea ) 
2 3 



Mavav (ea ) 








1 




2 




3 




4 




2.9 (- 


1) 


4.0 (- 


-1) 


4.2 (- 


-1) 


3.8 (- 


-1) 


3.0 (- 


1) 


4.2 (- 


1) 


3.3 (- 


-1) 


9.1 (- 


-2) 


1.2 (- 


-1) 


2.3 (- 


1) 


4.2 (- 


1) 


4.2 (- 


-2) 


2.5 (- 


-1) 


2.9 (- 


-1) 


1.7 (- 


1) 


3.3 (- 


1) 


2.2 (- 


-1) 


2.7 (- 


-1) 


2.5 (- 


-2) 


1.8 (- 


1) 


2.2 (- 


1) 


3.3 (- 


-1) 


4.9 (- 


-2) 


2.4 (- 


-1) 


2.2 (- 


1) 



M A ' u 'A"v" (ea ) 
2 3 



1) 2.8 (-2) 
4.9 (-1) 



3.1 (-3) 
2.0 (-2) 
6.0 (-3) 



2.2 (-4) 
2.0 (-3) 

2.0 (-2) 

6.1 (-D 



1.2 (-5) 
2.0 (-4) 
2.0 (-3) 
2.0 (-2) 
3.4 (-1) 





4.3 (- 



Mava'V (ea ) 
2 3 



2) 4.1 (-2) 
3.3 (-2) 



2.6 (-3) 

5.7 (-2) 
2.4 (-2) 



1.6 (-4) 
4.4 (-3) 
4.0 (-2) 
1.4 (-2) 



8.1 (-6) 

3.4 (-4) 
3.6 (-3) 

4.5 (-2) 
3.8 (-3) 



J"\J' 


1 

2 
3 
4 
5 



1.7 (-11) 
6.1 (-12) 
2.6 (-12) 
5.4 (-13) 
5.1 (-14) 



1 

7.9 (-12) 

1.8 (-11) 

3.5 (-12) 

4.6 (-13) 

6.9 (-14) 



2 

3.0 (-12) 
2.0 (-11) 

1.3 (-11) 
1.5 (-12) 
8.3 (-14) 



(cm 3 s l ) 
3 

2.3 (-12) 
7.0 (-12) 
2.3 (-11) 

9.8 (-12) 
7.0 (-13) 



4 

1.3 (-12) 

2.4 (-12) 
7.2 (-12) 
2.6 (-11) 

6.8 (-12) 



5 

4.5 (-13) 
1.3 (-12) 

1.5 (-12) 
6.9 (-12) 

2.6 (-11) 





6.1 (-11) 

1.2 (-10) 
3.9 (-11) 

2.0 (-11) 

1.1 (-H) 



1 



Z[C]fc y c ,,„ (em's- 1 ) 
2 3 



2.1 (-11) 2.8 (-11) 7.8 (-12) 

4.1 (-11) 4.0 (-11) 

6.1 (-11) 4.6 (-11) 

7.2 (-11) 5.5 (-11) 

2.7 (-11) 5.3 (-11) 4.9 (-11) 

1.1 (-11) 2.0 (-11) 4.3 (-11) 



4 

3.8 (-12) 
1.5 (-11) 
4.3 (-11) 
4.7 (-11) 

4.5 (-11) 



5 

2.3 (-12) 

6.4 (-12) 

1.7 (-11) 

4.5 (-11) 

4.8 (-11) 



2 [D]k D (cm 3 s- 1 ) 
2.29 (-13) 



/ ^(v)I v dv (s- 1 ) 
9.80 (-10) 



Z[D]k D (cm 3 s- 1 ) 
2.28 (-17) 



/ %o*(y)Iydv (s- 1 ) 
2.60 (-10) 



./ 
O 
© 







1 



n,/n(HCl) 
2 3 



10 3 

I 10 s 

* 10 7 

> 10 2 

I io 3 

* 10 4 

> 10 2 



10 3 

10 4 



1.0 (+00) 4.9 (-04) 8.5 (-06) 8.1 (-07) 6.1 (-08) 2.7 (-09) 

1.0 (+00) 4.9 (-04) 8.6 (-06) 8.2 (-07) 6.3 (-08) 2.9 (-09) 

1.0 (+00) 4.9 (-04) 8.7 (-06) 8.2 (-07) 6.3 (-08) 2.9 (-09) 

1.0 (+00) 7.5 (-04) 2.3 (-05) 8.9 (-07) 6.3 (-08) 2.9 (-09) 

9.7 (-01) 2.6 (-02) 1.4 (-03) 1.3 (-05) 2.2 (-07) 4.3 (-09) 

1.0 (+00) 5.0 (-04) 9.0 (-06) 8.5 (-07) 6.4 (-08) 3.0 (-09) 

1.0 (+00) 6.1 (-04) 1.3 (-05) 1.0 (-06) 7.7 (-08) 3.4 (-09) 

1.0 (+00) 1.7 (-03) 4.8 (-05) 3.1 (-06) 2.1 (-07) 8.4 (-09) 

1.0 (+00) 6.5 (-04) 1.4 (-05) 1.5 (-06) 1.5 (-07) 1.4 (-08) 

1.0 (+00) 2.1 (-03) 6.3 (-05) 7.4 (-06) 9.4 (-07) 1.1 (-07) 

9.8 (-01) 1.6 (-02) 5.6 (-04) 6.7 (-05) 8.9 (-06) 1.1 (-06) 



»//n(CO) 

1 2 3 4 5 

7.3 (-02) 2.5 (-01) 3.3 (-01) 2.2 (-01) 8.9 (-02) 2.9 (-02) 

7.3 (-02) 2.5 (-01) 3.3 (-01) 2.2 (-01) 8.9 (-02) 2.9 (-02) 

7.3 (-02) 2.5 (-01) 3.3 (-01) 2.2 (-01) 8.9 (-02) 2.9 (-02) 

7.1 (-02) 2.4 (-01) 3.3 (-01) 2.2 (-01) 9.4 (-02) 3.1 (-02) 
3.3 (-02) 9.8 (-02) 1.6 (-01) 2.0 (-01) 2.0 (-01) 1.5 (-01) 

7.3 (-02) 2.5 (-01) 3.3 (-01) 2.2 (-01) 9.0 (-02) 3.0 (-02) 

7.0 (-02) 2.3 (-01) 3.3 (-01) 2.3 (-01) 9.8 (-02) 3.3 (-02) 

5.3 (-02) 1.7 (-01) 2.6 (-01) 2.6 (-01) 1.6 (-01) 7.0 (-02) 

7.2 (-02) 2.4 (-01) 3.3 (-01) 2.2 (-01) 9.2 (-02) 3.1 (-02) 

6.3 (-02) 2.1 (-01) 3.1 (-01) 2.3 (-01) 1.1 (-01) 4.2 (-02) 
3.7 (-02) 1.2 (-01) 2.0 (-01) 2.4 (-01) 1.9 (-01) 1.2 (-01) 



( 1200 ll) : lFlowerl ( 120011) 



Rank etal] (fl965l); iTilford et all ( Tl970h ; Ivan Dishoeck et alj <fl982h; lLietal.1 JMTlb: iNeufeld & Greenl fl^h^data from 
~~ Eidelsberg et alJ dl992h : IChan tranupong jl992h : ISundholm et al.1 jl995l) : lGoorvitchl dl994h !rkirivama & Raol 



" data from 

iMacris &"P etropoulos J1985I 
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